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In two-dimensional systems possessing a high degree of symmetry, the repulsive 
electron-electron interaction produces a pairing force; the mechanism would fail in the 
presence of strong distortions. We have studied this in the one-band and three-band 
Hubbard Model. From partially occupied orbitals one obtains pair eigenstates of the 
Hamiltonian with no on-site repulsion (the W — pairs). The concept of W = pairs 
allows to make qualitative and quantitative predictions about the behaviour of inter- 
acting many-body systems, a quite remarkable and unusual situation. Exact numerical 
solutions for clusters with "magic" hole numbers reveal attraction between the holes 
in W = pairs. The effect occurs in all fully symmetric clusters which are centered 
on a Cu site; then holes get paired in a wide, physically relevant parameter range and 
show superconducting quantization of the magnetic flux. A canonical transformation 
of the Hamiltonian, valid for clusters and for the full plane, leads to a Cooper-like 
equation for the W — pairs. We have evaluated the effective interaction and found 
that W = pairs are the bare quasiparticles which, once dressed, become two-hole 
bound states. We applied the above theory to the doped antiferromagnet, and found 
that the ground state at half filling is the singlet component of a determinantal state. 
We write down this determinant and the ground state wave function explicitly in terms 
of a many-body W = eigenstate. Our analytical results for the 4x4 square lattice at 
half filling and with doped holes, compared to available numerical data, demonstrate 
that the method, besides providing intuitive grasp on pairing, also has quantitative 
predictive power. 



I. INTRODUCTION 

Most superconducting oxides [1] have large unit cells containing four or more elements and Cu02 planes 
separated by intermediate layers which stabilize the lattice and act as a charge reservoir. A large amount of 
experimental and theoretical work has clarified that mobile holes on these planes are responsible for super- 
conductivity. In the undoped solid the electrons are so strongly correlated that they cannot move freely as in 
a metal (antiferromagnetic Mott insulator). When holes are introduced into the plane by doping, the system 
becomes conducting and even superconducting provided the hole density is within the appropriate range. 
Increasing experimental evidence obtained in several cuprate superconductors has convincingly demonstrated 
that the pairing state of these materials has a non-trivial symmetry, ranging from s- to d-wave [2]. Moreover, 
the pairs exist above the critical temperature either in the form of superconducting fluctuations or preformed 
pairs. The latter aspect is apparent in the underdoped region in which a clear pseudogap essentially of the 
same magnitude as the superconducting gap is measured in ARPES [3], NMR [4], neutron scattering [5] and 
specific heat measurement [6] in metallic phase. All these signatures put strict constraints to any microscopic 
model of the cuprates. Any theory of the paired state must predict the correct symmetry and doping depen- 
dence of the binding energy of the pair, but the pairing mechanism must be doping-independent and robust 
enough to survive superconducting fluctuations well into the normal state, far from optimum doping. 

The class of theories based on correlation mechanisms start from single or multiple bands Hubbard-like 
Hamiltonians containing nearest neighbors hopping matrix elements t between Cu and O and 0-0 and large 
on-site Coulomb repulsive energies U between holes of opposite spins which tends to reduce double occupancy 
of these sites. For U » t, such Hamiltonians map into Heisenberg-like Hamiltonians for holes propagating 
through an antiferromagnetic background with an exchange interaction J between neighboring spins. Pertur- 
bative approaches such as the t- J Model [7] [8] , by means of the Guzwiller projection [9] force the holes to avoid 
on site interactions by letting U to be very large (or infinite) so that eventually the pairing may be achieved 
through residual small attractive interactions. In this way the problem is circumvented formally, but not 
understood physically. Such models do not identify from the outset the physical mechanism capable to nullify 
the large repulsive Coulomb interactions among the paired holes. On one hand, the Guzwiller projection cuts 



1 



a large part of the Hilbert space, and that would cost a catastrophic amount in kinetic energy, compared to the 
energy gain involved in the normal to superconductor transition, unless the repulsion is orders of magnitude 
larger than the hopping, which is not the case; on the other hand, the same idea could equally well be applied 
to a potato, i.e. it does not explain what is special in the CuC>2 plane. 

Due to the small size of the Cooper pairs, we believe that any serious proposal of a pairing mechanism should 
essentially be size-independent, i.e. it should work in the thermodynamic limit but it should be manifest in 
tiny systems as well. Exact calculations on finite models should bring to light interesting local aspects of 
the microscopic pairing mechanism, and be useful as tests for the analytic developments; pairing should be 
possible using a physically plausible range of parameter values. 

In this review article, we illustrate the "W = 0" pairing mechanism, which is specific of the square 2d 
structure of the CuC>2 plane; instead of forcing the system by the Guzwiller projection we get the attraction 
by using the symmetry. The symmetry provides a physical mechanism capable of cutting down the on-site 
repulsion from the outset. We want to show that the holes avoid the on site interaction via a correlation 
mechanism which is effective for any U/t. We demonstrate this correlation effect, showing that its basic 
ingredient is a configuration mixing driven by the symmetry of the system which entails the presence of 
degenerate hole states at the Fermi level. Moreover, we provide evidence that the same configuration mixing, 
when applied to the many-body problem, also produces pairing between the holes and superconducting flux 
quantization when physical parameter values are used. 

The work is divided into two main parts. In the first part, that is Section II, we review the pairing mechanism 
in small symmetric clusters. We begin by showing that by symmetry the Hubbard Hamiltonian admits two- 
hole singlet eigenstates that do not feel the on-site Hubbard repulsion, that we call W = pairs. First, we 
demonstrate that the prototype CuCU cluster already contains the essential ingredients of the mechanism. 
Then we prove a general theorem on the allowed symmetries of W = pairs. From the theorem we extract 
a practical technique to build them explicitly for symmetric clusters and for the full plane. We will see that 
no W — pairs exist if the symmetry is too low. In Section II B we show that the peculiar properties of 
W = pairs cause a stabilization of the interacting system which leads to pairing. A careful analysis is made 
on the smallest "allowed" cluster, the CuCU; then, the investigation is extended to larger fully symmetric 
clusters. It is shown that pairing can be obtained in a physical parameter range any time the number of 
holes is "magic" . Section II C is devoted to the analytic interpretation of the results. It is shown that in 
second order perturbation theory the effective interaction between the holes of a W = pair comes out from 
the interference of positive and negative contributions, depending on the point symmetry of the electron-hole 
states exchanged, and that it may be attractive. Finally in Section II D we demonstrate the supercounducting 
flux quantization and conclude that W — pairs are the bare quasiparticles which become two-hole bound 
states when dressed by the interaction. 

In the second part, that is Section III, we generalize the theory to the full plane. In Section III A, we review 
the band structure of the three-band Hubbard Model and fix some notations. In Section IIIB we show how 
W — pairs in the full plane arise at the Fermi level for any filling, and that their symmetry properties naturally 
lead to superconducting flux-quantisation. Section III C is devoted to the general canonical transformation 
leading to the effective Hamiltonian for pairs. The method is new and free from the limitations of perturbation 
theory; the relation of the present formalism to Cooper theory from one side and to the cluster results from 
the other is discussed. The final result of the analysis is an integral equation for pairs in the full plane, which 
is suitable for a numerical study. Therefore, in Section HID the integral equation is solved numerically. We 
first solve it in finite supercells and then take the asymptotic limit. Two kinds of bound states of different 
symmetries result, and the dependence of A on the filling and other parameters is explored. Next, we study 
the Hubbard Model at and near half filling. Here the problem is that the non-interacting ground state at half 
filling is highly degenerate, while the canonical transformation needs a unique non-interacting ground state to 
work. Therefore, in Section III E we remove the degeneracy in first order perturbation theory by means of a 
suitable formalism, the local formalism, together with the help of a theorem by Lieb. The half filled ground 
state at weak coupling is explicitly written down. We show that it is the spin singlet projection of a peculiar 
determinantal state; the latter which exhibits the antiferromagnetic property: the translation by a lattice step 
is equivalent to a spin flip. In Section IIIF we use the 4x4 square lattice as a test case. The half filled 
antiferromagnetic ground state is doped with two holes and an effective interaction between them is derived 
by means of the canonical transformation. The analytical results agree well with the numerical ones and this 
shows the predictive power of the approach. Finally, the main conclusions are summarized in Section IV. 
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TABLE I.: Hole parameters (in eV) for the three-band Hubbard Model in Eq.(l). 

II. PAIRING IN THE HUBBARD MODEL: SYMMETRIC CLUSTERS 

Our starting point is the three-band Hubbard Model 

H =K + V^Hubbard + W oH - sitc (1) 

where 

K = ^2 tas'ctaCa'a = tpd ^ {Cja ™ + U - C + l PP ^/ ( C J<7 C j'<7 + h.C.) + e d ^ h ia + € p ^ flj a (2) 
ss f ,a (ij) ,a (jj'),<7 

and 

WHubbard = Ud ^ niffin + ?7p njfhji, Woff-site — U pd ^ h i<y n ja , . (3) 

i i {ij),aa' 

Here, Cj (c^) are fermionic operators that destroy holes at the Oxygen (Copper) ions labelled j (i) and h — c'c 
is the number operator, (ij) refers to pairs of nearest neighbors i (Copper) and j (Oxygen) sites. The 
hopping terms corresponds to the hybridization between nearest neighbors Cu and O atoms, and are roughly 
proportional to the overlap between localized Wannier orbitals. Ud and U p are positive constants that represent 
the repulsion between holes when they are at the same Copper (d) and Oxygen (p) orbitals, respectively. U p d 
has a similar meaning, i.e. it corresponds to the Coulombic repulsion when two holes occupy two adiacent 
Cu-O sites. The on site energies e p and ed represent the difference in energy between the occupied orbitals of 
Oxygen and Copper. In the strong coupling limit, and with one particle per unit cell, this model reduces to 
the spin Heisenberg Model with a superexchange antiferromagnetic coupling [10] [11]. 

The Hamiltonian in Eq.(l) shows that for A p d = e v — ed > 0, the first hole added to the system will 
energetically prefer to occupy mostly the d-orbital of the Copper ions. This is indeed the observed situation 
in the undoped materials which have one hole per unit cell. When another hole is added to this unit cell, and 
working in the regime where Ud is larger than A p d, the new hole mainly occupy Oxygen orbitals [12]. This is 
in agreement with Electron Energy Loss Spectroscopy experiments [13]. From a band structure and by best 
fitting the results of up-to-date ab initio calculations [14] one can roughly estimate the actual values of the 
parameters in the Hamiltonian of Eq.(l). Our preferred set 1 is shown in Table I (in eV). 

Several years ago, Hirsch and Scalapino [15] introduced the quantity A(N) = E a (N)+Ef_)(N—2)—2Eo(N—l), 
where Eo(N) is the ground state energy of the system with N holes. They showed by direct diagonalization of 
several clusters that A can be negative with all positive Coulomb repulsion parameter. The A < behavior 
was taken as an indication that some kind of attraction was developing from repulsive interaction. From those 
calculations [16], it was apparent that negative A arose from unphysically large Cu-0 off-site interactions U p d 
of the order of 4-6 eV at least. Moreover the A < criterion has also been questioned [17] because when one 
introduces the nuclear degrees of freedom, the sign of A can be reversed. This can happen when the ground 
state with N — 1 holes is degenerate, and Jahn- Teller distorts, gaining energy by the distortion; the pairing 
then looks like an artifact due to the neglect of vibrations. Ultimately, the argument runs into trouble because 
it depends on a comparison of ground state energies with different hole numbers TV. 

In 1995, Cini and Balzarotti [18] resumed the subject from a new standpoint which avoids the need for 
a large U p d and is free from the above criticisms. They started from the observation that highly-symmetric 
clusters possess 2-holes singlet eigenstates which do not feel the on-site repulsion WHubbard, W = pairs. This 
property poses severe limitations on finite clusters that not only must possess the full Ca v symmetry, but must 
be centered around a Cu site, see Figure 1. 



1 We are using a symmetric gauge; for instance, F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988) 
use a different sign prescription for t p d and t pp , which corresponds to changing the sign of some orbitals; using 
their convention, the orbital symmetry labels ai and 6i get exchanged. 
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FIG. 1.: C11O4, C115O4 and CU5O16 clusters, the "allowed" geometries. Full dots are Cu ions. 



Such W = eigenstates are generally excited states of the two-body system and arise when the two holes 
occupy degenerate states of E = (x,y) symmetry. The interesting situation arises when the hole population is 
such that the non interacting ground state would have two E holes. Thus, any allowed cluster has its "magic" 
number of holes. The interacting ground state must be described in terms of a many-body configuration 
mixing, but the presence of the W = pair imparts to the many-body state a special character, as we shall 
show; the persistence of the "W = 0" character in the interacting case is the main point that we want to make. 
The need for a partially filled shell restricts the hole numbers for small clusters, whereas of course any even 
number of holes works in the thermodynamic regime. In the infinite plane, W = pairs are easily obtained 
[19] from the non-bonding band; a method for building them at any kinetic energy is discussed below. 

Consider the process of filling the one-body levels of the noninteracting system according to the aufbau 
principle. Since the order of the one-body energy levels depends on the site energies and hopping integrals, 
for a given cluster the "last pair" is a W = pair in part of the parameter space. Therefore, W = pairs fix 
the cluster structure, occupation and the range of the parameters. These restrinctions are so stringent that 
these clusters with such properties had not been studied previously. In particular, the present discussion does 
not apply to the forbidden geometries, like those examined by Hirsch et al. [15] [20] and Balseiro et al. [16]. 
The CU4O4 geometry considered by Ogata and Shiba [21] has the C\ v symmetry of the lattice, but lacks the 
central Cu, and therefore it is forbidden (the pairs on degenerate levels feel the on-site repulsion, as one can 
see by performing the two-hole configuration interaction calculation [18]). The dynamics of holes in CU4O4 
does not lead to pairing if the on-site repulsion on Oxygen is included [22] . 



In this Section, we limit the discussion to on-site interactions (U p d = 0), while the effect of the off-site terms 
is deferred to Section II B. 



For the sake of simplicity, let us first consider the CUO4 planar cluster of Figure 1 with 4 holes, one intrinsic 
and three due to doping. Site 1 is Cu, and the others are Oxygens. The single-hole energy level scheme vs. t pp 
is shown in Figure 2. The levels are labelled according to the irreducible representations of the dv symmetry 
Group whose character table is shown in Table II. 

At tpp = —0.65 eV, the first two holes go into a bonding level of A\ symmetry, and the next two go into a 
non-bonding level of E symmetry, with orbitals transforming like the (x,y) vector. Denoting the bonding state 
by a and the degenerate non-bonding states by e x and e y , one could build the degenerate singlet configurations 
|at a l e xT e xJ.) an d \ a 1 a l e y1 e yl) an d the triplet components |of 0|e x f e y f ) and \a-faie x ie y i) . Since parallel spin 
holes are not interacting in the model, in the Hartree-Fock (HF) approximation the triplet would have lower 
energy according to the Hund's rule. 

Yet, the singlet energy can be much improved by configuration interaction (CI). First, let us replace the 
(e x ,e y ) pair by the (a, f3), where a = (e x + e y )/\/2 and [3 — (e x — e y )/\/2 are oriented in the 0-0 directions 
and are connected by a 7r/2 rotation. The singlet configurations then read |*c«a) = \a^aict-\cn) and \^ pp) = 
\a^aiPif}{), both with HF energy 



A. W = Pairs 



1. W = Pairs in the CuOa Planar Cluster 




(4) 




FIG. 2.: One-hole energy levels of the CUO4 cluster vs. t pp ; the other parameter are specified in Table I. The 
levels e and fei are non-bonding, the a\ levels form a bonding-antibonding pair. The levels are labelled according to the 
irreducible representations of the Cuv symmetry Group. 
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TABLE II.: Character table of the C4 V symmetry Group. Here 1 denotes the identity, C 2 the 180 degrees rotation, 
C4 the counterclockwise and clockwise 90 degrees rotation, a x , a y the reflection with respect to the y = and 
x = axis and tr+, a— the reflection with respect to the x = y and x = — y diagonals. The C± v symmetry Group has 4 
one-dimensional irreducible representations (irreps) Ax, A2, B\, B 2 and 1 two-dimensional irrep E. In the last column 
it is shown how each of them transforms under C^, v . 



where s denotes the generic site and U s = Ud if s is a Cu site and U s = U p otherwise. The first term on the 
right hand side contains the one-body contributions 

K V a =)^(pl(s)t sa >(p a (s'), lf(s) = (s\(p) , If = a, Q (5) 

where t ss i is the hopping matrix element of the three-band Hubbard Model with diagonal entries t sa = 
(e p ) if s is a Cu (O) site, and n a (s) = \a a (s)\ 2 + \a a (s)\ 2 = n a (s) + n a (s). The configurations are mixed by 
{^cea\H\^ pp) — U s n a (s)nf3(s). The lowest eigenvalue corresponding to the choice 

= -L(l* aa >_l* /s/J » (6) 

is 

e = (y aa \H\y aa ) - (*„a|H|* w > = ^(K„ + Kccc) + J2 U *i n *( s ) + 2n a (s)n Q (s)] + W, (7) 

a s 

where W = ^ U 3 n a (s)[n a (s) — np(s)] represents the direct interaction between the non-bonding holes. At 
this point the issue is decided by the peculiar property of the system that n a (s) and np(s) are identical, despite 
the fact that the quantum states are distinct and orthogonal, and therefore the interaction between those holes 
is turned off exactly. We emphasized that the existence of degenerate states such that n a (s) = np(s) is the 
only condition. This two-hole singlet eigenstate with vanishing on site repulsion is a W = pair. The above 
argument brings out the essential role of symmetry, and can be extended to clusters of arbitrary size. The 
singlet in Eq.(6) turns out to be 1 B 2 which transforms as xy and has a totally symmetric charge distribution. 

Since W = for both 1 B 2 and the triplet, two-body calculations cannot predict which is the ground state. In 
Section II B we shall show that pairing is indeed achieved in part of the parameter space and a non-degenerate 
singlet ground state prevails. 
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2. W = Theorem 



For any finite cluster having the full C± v symmetry and a Cu ion at the center, one can obtain W — pair 
eigenstates [23] even if the above procedure is not directly extendible. 

Here we want to point out a more powerful and elegant criterion to get all the W = pair eigenstates 
of a given system. As we shall see this criterion applies to Hubbard Models defined on an arbitrary graph 
A. Let Q be the the symmetry Group of the non-interacting Hubbard Hamiltonian K — Y] x gA a t xv c xa c va . 

By definition, every one-body eigenstate of ^Hubbard = K + Wkubbard, Wkubbard = J2xeA U x n x ^h xl , can be 
classified as belonging to one of the irreducible representations (irreps) of Qo- We may say that an irrep r\ is 
represented in the one-body spectrum of i?Hubbard if at least one of the one-body levels belongs to r\. Let £ be 
the set of the irreps of Qo which are represented in the one-body spectrum of -ffHubbard- Let \ip) be a two-body 
eigenstate of the non-interacting Hamiltonian with spin S z = and P^ the projection operator on the irrep 
r). We wish to prove the 

W = Theorem. - Any nonvanishing projection of \ip) on an irrep not contained in £ , is an eigenstate of 
-f^Hubbard with no double occupancy. The singlet component of this state is a W = pair. Conversely, any 
pair belonging to an irrep represented in the spectrum must have non-vanishing ^Hubbard expectation value: 

■n^£^Wuu bb , ld P M W=0 (8) 

Proof: Let us consider a two-body state of opposite spins belonging to the irrep r/ of Qo- 

l^) >= £v (, W)4 T c lil°>- ( 9 ) 

xy£A 

Then we have 

h x1 h xl W M )=iP M (x,x)cl 1 ci i \0) = ^\x,x)\x],xi). (10) 
We define P^ as the projection operator on the irrep r\. Since 

P MJ2^ M (x,x)\x T,x |) = Y,i> M (x,x)\x \,x I) (11) 

if P^\x |, x I) = Va; G A then xp^ (x,x) = Vi £ A. It is worth to note that this condition is satisfied if 
and only if P M \xcr) = Vx G A where \xa) = cJ. CT |0). 

Now it is always possible to write \xa) in the form \xa) = 6£ {x)\r) a ) where \r) a ) is the one-body 

eigenstate of //nubbard with spin a belonging to the irrep r\. Hence, if rj ^ £ , P^ ^\xa) = and so P«)\ x t 
, x I) = 0. Therefore, if |V>^) is a two-hole singlet eigenstate of the kinetic term and r\ £ £ , then it is also an 
eigenstate of W-Hubbard with vanishing eigenvalue, that means a W = pair. 

There are cases when this theorem puts useful restrictions on the possible ground state symmetries; this 
happens when the ground state with N holes corresponds to filled shells in the U — *■ limit and we wish to 
study the N + 2 hole ground state. At least for small U the two added holes will form a W = pair, since 
this minimizes the ground state energy. Let Q C Qo be the symmetry Group of // Hu bbard and U C (N) the 
minimum crossover value of U, that is, the ground state with N holes \$u(N)) has well defined symmetry for 
< U < U C (N). Let \& u= o{N)) be non-degenerate (e.g. closed shells case) and belong to the one-dimensional 
irrep rjo(N) G Qo- When we add the two extra holes, the new ground state for U = is a, W — pair over 
\®u=o{N)) , and its symmetry is r)o(N) ■ r)w=o, where rjw=o G Qo is the symmetry of the added pair. Turning 
on the interaction r)o(N) ■ r\w=o breaks into the direct sum of irreps belonging to Q and among them there 
must be the symmetry of \<&u(N + 2)) if < U < U C {N + 2). This remarkable restriction is posed by Group 
Theory alone: which of the symmetries that remain allowed is actually realised in the ground state depends 
on the dynamics. 

The complete characterization of the symmetry of W = pairs requires the knowledge of Qo- A partial 
use of the theorem is possible if one does not know Qo but knows a Subgroup Qq ■ It is then still granted 
that any pair belonging to an irrep of Qq not represented in the spectrum has the W — property. On the 
other hand, accidental degeneracies occur with a Subgroup of Qo, and by mixing degenerate pairs belonging 
to irreps represented in the spectrum one can find W = pairs also there. 

The theorem tells us that the bigger is Qo, the larger is the number of W = pairs. Indeed, for a given 
graph, the number of one-body eigenvectors is fixed, while the number and the dimension of the Qo irreps 
grow with the order of the group. Therefore, also the number of irreps not represented in £ grows, and this 
means more W = pairs. At the same time, a big "non-interacting" symmetry Group Qo grants large level 
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-1.65603 
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-2.88605 




2.21855 


So (4) 


8.52 


7.01287 
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1.46 


-0.064 


0.00747 



TABLE III.: Ground state energy Eq(N) of the C11O4 cluster with N = 1 to 4 holes and effective repulsion A (4) = 
Eq(4) — 2Eq(3) + Eo(2) between the holes of E symmetry in eV. HF: Hartrce-Fock approximation; CI: configuration 
interaction; Exact: numerical diagonalization of the Hamiltonian. The parameters are those reported in the Table I 
except for U p d that is taken to be zero. 

degeneracies. As we have seen, this fact allows the existence of W = pairs formed by degenerate orbitals. 
In the next Section we show that an anomalously low effective repulsion takes place in the interacting system 
when such levels are on the Fermi surface and that in a certain region of the parameters space this leads to 
pairing. 

Finally we observe that the W = theorem predicts W = pairs of A2 and/or B2 symmetry in the CuCU 
cluster, as one can see by inspection of Figure 2. This is consistent with the above results where the symmetry 
of the W — pair in Eq.(6) was estimated to be B-z- Further, when t vv = 0, the non-bonding levels e and 61 
in Figure 2 become degenerate and a new W = pair of Ai symmetry is obtained. At a first sight, this fact 
seems to contradict the W = theorem, but it is not so. As we have observed before, unexpected symmetries 
arise if one use a Subgroup of Go- For t pp — 0, the extra degeneracy cannot be explained in terms of the C 4v 
Group, whose irreps have dimensions 2 at most. It is not hard to realize that for t pp — any permutation of 
the four Oxygen sites is actually a symmetry and therefore the symmetry Group is enlarged to II4 (the group 
of permutations of four objects) which admits irreps of dimension three. Then, labelling the one-body levels 
with the irrep of II4, one readly realize that irreps not contained in the spectrum give W = pairs, according 
to the W = theorem, and that among them there is one containing the A\ irrep of Cav 

B. CI and Exact Solutions for Symmetric Clusters 

In the present Section, we wish to show how A < may be induced by symmetry through degeneracy. 

1. A < in the Cu0 4 

Let us first consider the CUO4 cluster with four holes in the non-interacting limit: according to the aufbau 
the first two holes go into a bonding level of Ai symmetry and the next two into a non-bonding level of E 
symmetry. Group Theory predicts that the interactions resolves the sixfold-degenerate ground state E <g> E 
into 3 ' A2, 1 Ai, 1 B\ and 1 £>2- The two- body calculations say that 1 Bi and the triplet are lowest, since W = 
for both. The ground state will turn out to be singlet, while Hund's rule would have predicted 3 A2; we shall 
see why and how low spin prevails. 

In Table III we report the ground state energies Eo(N) of the CUO4 cluster with N holes, N = 1 to 4, and 
the effective repulsion energy A(4) = Eq(4) — 2Eq(3) + Eo(2). While W gives the direct interaction between 
the two holes in the CI calculation, A(4) contains indirect interaction as well. The HF results correspond 
to a spontaneously broken symmetry of the 4-hole singlet and fail to reproduce the quantum effect we are 
discussing. Accordingly, the effective on-site repulsion is A(4) = 1.46 eV, which is quite a large repulsive 
barrier to overcome before pairing can be attained. In the CI calculation it is assumed the form (6) for the 
wave function and optimized the energy with respect to the orbitals. This does lead to a small negative 
estimate for A (4). The exact diagonalization with the same parameters yields a positive but tiny A (4) with a 
1 i?2 ground state (the symmetry of the W = pair), while the first excited state is a 3 A2 triplet lying D — 18 
meV above it. This suggests the interpretation that W = two-hole eigenstates are bare quasiparticles and 
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-0.0049 
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-0.0364 
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TABLE IV.: Effective interaction A(4) and energy D of the first excited triplet relative to the singlet ground state 
as a function of t pp for the Cu04 cluster. The off-site U p d parameter is still ignored. 

their attractive screened interaction produces the pairing 2 . Although the CI calculation gives a negative A (4) 
value, the exact calculation yelds a tiny positive one. However, if t pp is increased, A(4) does become negative, 
which shows that the correlation effects may lead to pairing (see Table IV). Even if the estimation of the 
singlet-triplet gap depends on the parameters, it is consistent with the 45 meV excitation observed in High 
Resolution Electron Energy Loss spectra [24] and the 41 meV resonance found in Polarized Neutron Scattering 
experiments [25] on superconducting YBa2Cu307. 

The maximum binding energy occurs at t pp = 0, when all the non-bonding orbitals are degenerate and 
the configuration mixing is thereby enhanced. On the other hand, at positive t pp 's, A(4) becomes large and 
positive, because the bi non-bonding level is pushed below the degenerate one, and the mechanism is hampered 
(at tpp = +0.65, A(4) = 0.53 eV). 

Next, we discuss the dependence of A (4) on the other parameters [18]. If we decrease e p , A(4) decreases 
because this makes the system more polarizable. The e p dependence when all the other parameter are kept 
fixed and U p d = is almost linear down to e p — 0. By contrast, the dependence of A (4) on the Cu-0 hopping 
parameter t p d is complex. For small t p d, A(4) is large and positive, then falls to zero at t p d ~ 0.7 eV. Negative 
A(4) values are obtained for t p a > 1.4 eV. Finally let us consider the effects of a non-vanishing U p d- According 
to Ref. [16] positive U p d values do not spoil the mechanism, and tend to be synergic with it. Indeed, values 
of Upd > 0.6 eV give negative A(4) values even for t pp = —0.65 eV (in the range 0.2-1.2 eV considered in Ref. 
[18] A (4) is a monotonically decreasing function of U p d). 

As we have seen, symmetry plays a key role in the analytical calculations; the numerical results on the 
CuCU cluster further demonstrates the fact that any lowering of the symmetry is reflected by a corresponding 
increase of A(4) [18]. This implies that the CI mechanism is working thanks to the C\ v symmetry, and it is 
cut off whenever the symmetry is lowered. 

2. A < in Symmetric Clusters 

The total hole concentration which corresponds experimentally to the superconducting state is ph ~ 0.4 
per atom. In the CUO4 cluster we need four holes to operate the CI mechanism, because we need holes in 
degenerate states. Therefore p h is too large by a factor of two, and the W = pair involves nonbonding states, 
which are far from the Fermi level in the bulk. It is important to realize that these undesirable features are 
peculiar of the prototype CUO4 cluster, and already disappear in CU5O4, the next larger cluster of the same 
symmetry, see Figure 1. In fact, four holes are still sufficient to reach degenerate states, as shown in Figure 3 
(a), but ph ~ 0.44 is much closer to the experimental value. 

The degenerate states are bonding, in contrast with the CUO4 case, and the holes now spend a large fraction 
of their time on Cu sites. The effective interaction A(4) for CU5O4 and the spin triplet separation D vs. t pp is 
shown in Figure 3 (b). A (4) is negative in a wide parameters range by a correlation mechanism that does not 
need off-site U p d to operate. To visualize the physical effects of the effective attraction, one can compare the 
ground-state charge distribution for the non- interacting and interacting holes [23]. When A(4) < 0, the charge 
distribution shrinks; besides, the holes tend to accumulate more on the Oxygens where their mutual repulsion 



2 We note that the experimental evidence on the symmetry of the pairs is still conflicting, and different 
techniques tend to give different results. Angular resolved photoemission data favour d pairs with an order 
parameter of 1 Bi(x 2 — y 2 ) simmetry. The W = pair of symmetry 1 B2(xy) is consistent with these findings, 
since its total-symmetric charge density has nodes at 45 degrees from the Cu-O axis. Phases are less significant, 
and we can switch irreps by a gauge transformation. In the full plane the present approach also predicts W — 
pairs of different symmetries, and the symmetry of the most bound pairs depends on the filling. It is fair to 
say that this problem is far from being settled. 
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is mitigated by the CI mechanism. The opposite takes place in the "unpaired" region where the charge of the 
two least bound holes piles up preferentially on the external Cu sites. The size of the problem with four holes 
is 1296; it reduces to 81 on projecting to the B2 symmetry of the ground state, while the symmetry of the 
first excited state is a 3 ' A2 triplet. As in the Cu04 cluster the maximum binding energy occurs at t vv = 
[19] and takes the value A (4) = —15.7 meV; thus the binding energy |A(4)| is about half the CUO4 value. 
For larger clusters the computational effort is more demanding (44.100 configurations for the CusOie). For 
t pp = 0, and U p d = 0, it was found [23] that A(4) = —7.5 meV in the CU5O16 planar cluster (with U p d = 1 
eV, A(4) = —7.6 meV). The calculated trend of A(4) with the number of atoms |A| is shown in Figure 4 (a). 
The decline of ph with increasing number of atoms is slow enough to allow a significant contribution to the 

o 

binding energy for cluster sizes that correspond to the actual size (f ~ 10 — 20 A) of the superconducting pair. 




FIG. 3.: (a) One-electron energy levels in the CU5O4 cluster versus t pp . The levels are labelled according to the 
irreducible representations of the Cav Group. We find a couple of E symmetry (twofold degenerate) levels (one bonding 
and one antibonding), pairs of bonding-antibonding levels of A\ and B\ symmetry and a nonbonding level a\. (b) 
Dependence of A(4) and D on t pp for the CU5O4 cluster. In both (a) and (b) the parameters (in eV) are e p — td = 3.5, 
tpd = 1-3, Up = 6, Ud = 5.3 and U p d = 0. 

As shown in Figure 3 (a), there are two degenerate levels available to test the W = pairing mechanism; 
moreover electron pairs and hole pairs are related by an approximate charge conjugation symmetry and we 
can perform the test by filling the levels from bottom to top and the other way. The CU5O4 Hamiltonian was 
diagonalized [26] and the ground state energy was obtained as a function of the number N of holes. Eventually, 
we find electron pairing with a binding energy A(N) of a few tenth of meV in the physical parameter space. 

Let us consider the three-band Hubbard Hamiltonian of Eq.(l). By a canonical transformation from holes 
to electrons, c sa — > c| CT , we get for U p d = 

H = ^2,{2e s + U s )-^2 l {e 3 + U s )c\ a c sa - ^ t M >c\ a c a , a + ^ U s h s fh sl . (12) 

3 SCT (s,s')<7 S 

Thus, in the electron representation the signs of the site and hopping integrals are reversed. With respect 
to the hole case, the sequence of single-electron energy levels of the CU5O4 cluster is such that bonding and 
antibonding states are interchanged. In the hole representation the antibonding state of e symmetry contains 
two holes when 12 holes are present in the cluster. Therefore Nh~4, 12 are magic numbers (to avoid confusion 
we will denote by Nh {N e ) the number of holes (electrons)). In the electron representation W = electron 
pairs are obtained for N e — 4, 12. Letting (N e ) be the ground state energy of the cluster with N e electrons, 
the effective interaction between the two electrons in the pair is measured by 

A e (iV e ) = 4 e) (Ay + E^(N e - 2) - 2E^(N e - 1). (13) 

Since the dimensionality of the one-body basis is 18 

A e (4) = £< e) (4) + £< e) (2) - 2£< e) (3) = E ( h) (U) + E ( ' l) (16) - 2£<' l) (15) = A h (16) (14) 

where EQ h \Nh) is the ground state energy of the cluster with Nh holes. Similarly one finds A e (12) = Ah (8). 

We recall that we speak of electron pairs when two added electrons partially occupy a degenerate state and 
of hole pairs when the same situation is reached by adding two holes; however the final situation is exactly 
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100 




FIG. 4.: (a) A(4) and the singlet-triplet separation D for allowed clusters containing |A| atoms. The dots are the 
values calculated with the parameters listed in the Table I but t pp = 0. The curves are a guide to the eye. The clusters 
are Cu0 4 , Cu 5 4 , Cu 5 8 , Cu 5 Oi 2 and Cu 5 Oi 6 . (6) Pairing energy A(N) = E (N) + E (N - 2) + 2E (N - 1) of the 
CU5O4 cluster with N = 4 and N = 14 as a function of t pp . (c) A(jV) as a function of t pp for N = 10 and N = 12. 



N h 


4 


6 


8 


12 


14 


16 


A (meV) 


-15.7 


1469.2 


-10.85 


43.72 


1109.2 


-25.47 



TABLE V.: Exact diagonalisation results for A(N h ) (meV), using U p =6 eV and U d =5.3 eV. For TV,, =4, 8 and 16 
pairing takes place, and at 7V^=12 the repulsion is drastically reduced. For N^=6 and 14 the W = pairs are not 
involved and the normal repulsion develops. 

the same. One readily verifies that in a canonical transformation from holes to electrons the two-hole state 
becomes a two electron state of the same form. Therefore the two-body W = state is invariant under charge 
conjugation, and if holes are paired, electrons are also paired. In order to avoid switching all the time between 
the two equivalent pictures, below we discuss everything in terms of holes. 

In the symmetric CU5O4 cluster we previously found that the exact diagonalization of the Hamiltonian 
matrix with four holes in a wide range of negative t pp values for the 1 £>2 singlet gives a negative A (4). 
Therefore one expects that A(Nh) is small or negative for Nh = 4, 8, 12 and 16 and large and positive for any 
other Nh- 

The t pp dependence of A(Nh) computed using an enhanced Lanczos diagonalization routine [26] is shown 
in Figure 4 (6) and (c) for A(4) and A(16) and for A(10) and A(14), respectively. In Table V we summarize 
the results for A(Nh) when the model parameters are fixed according to the Table I but t pp = and for a 
wide range of Nh- One sees that for Nh ~ 4, 8, 12 and 16 A is much smaller in absolute value than for other 
fillings, as predicted. This confirms that W — pairs are involved. 

It is apparent that a considerable degree of symmetry under e-h exchange exists in the cluster and that 
the pairing interaction is not restricted to a single value of doping. The electron-hole symmetry is not exact, 
because the one-electron energy level spectrum of Figure 3 (a) is not invariant when up is exchanged with 
down and left with right (t pp — > — t pp and E — > —E); however one could have predicted the negative A and its 
magnitude by the approximate symmetry. In real systems electron superconductivity has been first reported 
[27] in the Nd2Cu04 doped with Ce with a maximum T c =24 K for optimal doping concentration 5=0.15, 
similarly to more common hole superconductors. Nd2_iCeiCu04 is an example of high-T c cuprate in which 
charged carriers are electrons. 

To summarize, in this Section we have shown that a correlation effect (namely, the CI mechanism) leads to 
negative A values by means of the exact removal of the direct Coulomb interaction between pairs of degenerate 
hole states. The CI mechanism is a consequence of the symmetry and electronic structure of the Cu0 2 plane, 
and do not depend on the size of the cluster. The concept of W — pairs is useful already at a qualitative 
level to make predictions about the behaviour of interacting many-body systems, and this is quite remarkable 
and unusual. Any large deviation from the full symmetry restores the normal on-site repulsion The binding 
energy A is of the correct order of magnitude (10 meV) for the set of parameters which corresponds to the 
best current estimates for the actual system. In the next Section we will definitively remove the ambiguity 
concerning the physical meaning of A. We will show that in fully simmetric clusters, as the ones studied 
above, negative A corresponds to pairing. 
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FIG. 5.: Arrangement of the two CUO4 planar cluster E and Ej, described in the text, t is the hopping integral 
between the two Cu sites. 



C. The Physical Meaning: Effective Interaction 

In this Section we show that in the Hubbard Model the negative-A behaviour due to the "W = 0" mechanism 
has the physical meaning of hole pairing, due to an effective attraction between the holes of the W = pair. 
The fact that pairing occurs can be estabilished by the following argument [23] . Consider a system E consisting 
of a couple of identical clusters (say, CuCU) that we denote E a and E b . Let E a and Ei, be placed one on top 
of the other and suppose the two systems are linked by a interplane hopping integral r << A(4) as shown in 
Figure 5. Such an arrangement is chosen arbitrarily for the sake of definiteness, but the important thing is 
that the two clusters are linked in such a way that the symmetry is not disturbed and the magnitude of r 
is such that its influence on the ground state energy is negligible compared to the other energies of the system. 
We may then write the full Hamiltonian as H — Hs a + Hs b + H T and regard the weak link term H T as a 
perturbation. In zeroth order, H can be diagonalized simultaneously with the occupation numbers Ns a and 
Ns b of the two clusters. The interesting situation arises when the number of holes in E is N = Ns a + Ns b = 6. 
The nature of the ground state depends on the parameters. The "normal" state is obtained if A(4) > 0; the 
unperturbed ground state is |5>normai) = \Ns a = 3} <8> \N^ b = 3), that is, the product of the individual ground 
state of two 2 E independent clusters containing three holes each. The unperturbed ground state energy is 
« 2£ (3). 

On the other hand, if A(4) < 0, that means Eq{4) + Eo(2) < 2_Bo(3), the zeroth order ground state is 

|$paired) = |Afc„ = 4) ®\N. b =2)±\N. a =2)*\N» b = 4) (l5) 

v2 

and the system prefers to disproportionate. It is clear that A(4) is just the energy needed to break the pair. 
The ± sign arises because the system has a bonding and an antibonding level; to lowest order, the two levels 
are separated by an energy of the order of r 2 /| A|. If the system is prepared with a pair of holes in one of the 
planes, it evolves in such a way that the pair will jump back and forth from E a to Ej, with that frequency 3 . 
This is clear evidence for pairing; if many systems like E are connected by weak links in a similar way, pairs 
will propagate 4 . 

Next, we want to show that in the Hubbard Model the pairing is due to an effective attraction between the 
holes of the W = pair. This point is more suble; one could argue that A < depends on a comparison of 
ground state energies of systems with different hole numbers N, and hence it is not obvious that it is a property 
of any of them. Moreover, the definition of A in terms of systems with different numbers of particles leads to 
a further difficulty. When one introduces the vibrational degrees of freedom, the state with a sigle hole in the 
degenerate level can gain energy by a Jahn- Teller distortion. The sign of A can even be reversed. The pairing 
then looks like an artifact due to the neglect of the vibrations. Therefore we want an alternative definition 
of A which is equivalent to the previous one at weak coupling but depends exclusively on the properties of 



3 We remark that interplane oscillations (Josephson plasma mode) had been observed in La^-iSraCuO-i, see 
Ref. [28]. 

4 Due to the mixing of degenerate states, the ground states of E belonging to positive and negative A(4) also 
belong to different symmetries. Indeed |$ pa irod) is nondegenerate, since \Ns a = 4) is 1 B2 while \N^ b = 2) is 
Mi. On the other hand, 2 E is fourfold degenerate and l^normai) has degeneracy 8 in the S z = subspace. 
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the Hubbard Model in the presence of the pair.' 1 We first show how to define the effective interaction in 
clusters, where we may take advantage from the exact diagnalizations. In all the allowed clusters up to 21 
atoms, the lowest one-hole level belongs to Ai symmetry, and the next e level yields the W = pair. In 
the previous Section we noted that the interactions produce a non-degenerate 1 i?2 4-hole ground state with 
A < 0. Below, it is shown [29] that in the case of dressed W — singlet pairs, A < arises from an effective 
pairing interaction which is attractive; the same interaction is repulsive for triplet pairs. The argument rests 
on a comparison of the perturbation series for A and for the two-hole amplitude; for the dressed W = pair 
the two-hole amplitude actually involves A and yields its dynamical interpretation. 



1 . Diagrams for A 

The perturbation series for A in powers of U is obtained by the well-known diagrammatic expansion of the 
ground state energy Eo of a many body system. This is based on the theorem 

E = £< 0) +i lim In R(t) (16) 

where E^ is the non interacting ground state energy, r\ = +0; 

R(t) = ($o|t/(t)|$ ), (17) 

|3>o) is the non interacting ground state, and U(t) is the time evolution operator in the interaction repre- 
sentation. One then expands the correlation function R(t) using the linked cluster theorem, that simplifies 
the expansion, but the diagrams that violate the Pauli principle and/or the number of particles must be 
retained. The diagrammatic rules are readily applied to the N — 2 case, when |$o) is a non degenerate 
single-determinant state. For N = 3 and N = 4, the non interacting ground state is degenerate, while the 
derivation of Eq.(16) assumes that the interactions do not modify the symmetry of the ground state; therefore 
we take |$ ) as the x component of the 2 E irreducible representation for N = 3 and 1 B 2 for N = 4; we know 
from direct diagonalization that these symmetries are correct for the whole series of clusters. For N = 4 the 
1 i?2 component of the non interacting ground state reads 

|$o) = !f£fg>+Jfggg> (18) 

where in terms of orbitals \e x e y ) is the state \a^aie x ^e yi ). Since this cannot be written as a single determinant, 
some care is necessary in the diagrammatic expansion; the correlation function is 

R(t) = {e x e y \U(t)\e x e y ) + {e y e x \U(t)\e x e y ) . (19) 

Applying the linked cluster theorem, the first term leads to the standard diagrams, while the second is 
anomalous, in the sense that it vanishes in the non interacting limit. Thus, one finds all the diagrams 
contributing to the expansion for U(t) averaged over |of aie x f e yi ), plus the anomalous or spin-flip diagrams 
that have entering e x f and e y i and outgoing e w f and e x | lines [29], see Figure 6 (A). 



5 This argument clarifies the nature of the pairing in the Hubbard Model, which is due to an effective attraction 
leading to a bound state. Including vibrational degrees of freedom in the description of a cluster or a plane, 
what the system would do with more or less holes than it actually has is unimportant if we can define the 
pairing by means of an effective interaction. Much more complex modeling would of course be necessary to 
study coupled subsystems like those of Figure 5 in the presence of vibrations. 
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FIG. 6.: (A) Linked diagrams contributing to the second-order expansion of Eq(A). Diagram on the left is normal, 
while the one on the right is anomalous. The a, b, e x and e y lines are labelled according to the irreps of the C*4„ Group. 
(B) Typical spin-flip diagrams for the anomalous propagator A s f. 

Let W(k,l,m,n) = (fc|/||WHubbard|n^T n l)i the anomalous diagram can be obtained from the normal one, 
which is proportional to W(a, b, e x , e x ) 2 , simply by letting W(a,b,e x ,e x ) — > W(a,b, e y , e y ) in the upper in- 
teraction line. A sign change follows, in agreement with the standard diagrammatic rules. In summary, one 
starts the expansion as if the ground state were the single determinant \a^aie x ^e Vl ); the diagrams that involve 
propagating e x and e y lines get corrections, which are obtained by exchanging e x and e y in the upper inter- 
action line where the e x or e y lines end; the loose lines should then be joined, and the diagram so obtained 
carries a minus sign because of the change in the number of loops. It is convenient to denote the hole orbitals 
by b for the 61 's and a (a 1 ) for the occupied (empty) orbitals 01. After some algebra, one obtains [29] the 
second-order approximation to A (4) 



A (2) (4) = -2 



\ - W(a,b,e x ,e x ) _ \ - W(a, a', e x ,e x ) 
4^ ( £ i> - e ") ( e a> ~ e a ) 



(20) 



where e a is the one-hole energy of the a orbital and so on; the sums run over all empty states of the appropriate 
symmetries (ai and 61), while no contributions arise from the empty e orbitals since the relevant IFmibbard 
matrix elements vanish. Eq.(20) holds for any symmetric cluster up to 21 atoms. The sign of A is seen to 
depend on the relative weight of the virtual transitions to states of different point symmetry, and ultimately 
on the parameters in the Hamiltonian. 



2. Diagrams for the Two-Hole Amplitude 

Let G2 denote the two-hole amplitude for holes of opposite spin in the degenerate (e x ,e y ) orbitals. Singlet 
and triplet arise from the space spanned by the non-interacting states \e x e y ) and \e y e x ), which are connected 
by Cav operators. In the Nambu formalism, G2 is a matrix 

G ^ = (l A a)-> (21) 

here A is a normal propagator, while A s f involves a spin-flip and vanishes in the non- interacting limit. When 
interactions are included, the first order contribution to the scattering amplitude vanishes because the product 
of e x and e y orbitals vanishes on all sites. Therefore, the spin-flip process, A s f, shown in Figure 6 (B), is 
the lowest-order e x — e y scattering, and produces the effective interaction to second-order. The system makes 
virtual transition to 4-body (3 holes and 1 electron) states. 

Evaluating A s t according to the standard rules, one finds [29] that near u> — 2e e , A s t w ,' a( ^ 4 A +less 

[uj 2t e ) 

singular terms; moreover, to second order, A = ^J 2e . We know from symmetry that G2(co) has singlet and 
triplet eigenvectors -^-(1,±1); for the singlet and triplet we get 

A + A,f = —!—, A-A sf = —!—, (22) 



where r) Sz t are the new eigenvalues. Therefore, for the singlet the expansion is 

1 1 A (2) (4) ^ 1 



u-ris LU-2e e (w - 2e e ) 2 w - 2e e - A( 2 ) (4) ' 



(23) 
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that is r/ s = 2e e + A' 2 ' (4); the triplet receives the opposite correction, and the 1 £>2 — 3 Ai separation is 
D = 2|A(4)|. Negative A(4) means that the spin-flip interaction is attractive for 1 B2 which is pushed down 
by |A(4)| and becomes the ground state, while the triplet is pushed up by the same amount. In this way, A(4) 
can be redefined without any reference to the ground state of clusters with a different numbers of holes, and 
we are free from the objections based on the Jahn- Teller distorsion of odd-iV clusters. 

Since U p and Ud are not small compared to the Cu-0 hopping term t p d, the second-order is generally a poor 
approximation; interestingly, W — pairs are an exception, because the large interactions are dynamically 
small. The first-order term vanishes, and the second order is of the order of tens of meV. Comparison with 
exact numerical diagonalization results shows that the second-order approximation for D and A(4) is already 
rewarding. In Figure 4 (a) it is shown how D and A(4) scale with the cluster size; they are indeed very 
closely proportional, although their ratio is somewhat larger than 2. Thus, the comparison with D is useful 
to show that A(4) really measures a property of the 4-hole state, and that the two E symmetry holes actually 
experience attraction. 



We test the paired state by exposing the system to a vector potential A(r) according to the Peierls pre- 
scription, 



where <j>o = 2-7r/e is the flux quantum, and looking for superconducting correlations. In a macroscopic exper- 
iment, one makes a sample with a large hole, and inserts a magnetic field; the ground state energy Eq(4>) is 
trivially a periodic function of 4>/4>o, where <f> is the flux in the tube, and <f> = is an extremum . Supercon- 
ductors quantize the flux by allowing integer and half- integer multiples of <fio, because the quasiparticles that 
screen the vector potential carry charge 2e. A flux 0o/2 corresponds to a superconducting ground state with 
pairs having different symmetry than those prevailing at (f> = 0, and the system is stable; arbitrary <f> values 
are not allowed because they cost large amounts of (free) energy. 

Canright and Girvin have shown [30] that cluster calculations can be used to obtain qualitative insight on the 
occurrence of superconductivity, by looking for a tendency to flux quantization. The signature [31] is present 
when <f> — is a minimum of Eo((j>) and the only other minimum of comparable depth occurs at (f)/(f>o = §; the 
barriers separating the minima are small in a small system but one expects them to increase with size, leaving 
the flux quantized in units of (f>o/2. Canright and Girvin [30] used a square lattice of rectangular geometry 
and periodic boundary conditions along one of the axes; to demonstrate the effect, they assumed an attractive 
on-site interaction and observed superconducting flux quantization for even hole numbers and strong enough 
attraction. In the present problem, with a repulsive Hubbard Model, the mechanism of attraction is driven 
by the Cav symmetry, and cannot operate with such an unsymmetric geometry. The flux must be inserted in 
such a way that the system is not distorted. We must insert the flux tube inside allowed clusters, but since no 
closed path encircling the central Cu is available for the holes, there is no response to a central flux tube. Cini 
and Balzarotti [29] used an external device providing a closed path around the flux tube. The device, however, 
cannot be chosen arbitrarily, because the cluster must remain allowed and the singlet pair must remain one 
with W — 0. To fulfill this condition, Cini and Balzarotti introduce an infinitesimal hopping tdd between the 
external Cu's (see the inset of Figure 7) and study the linear response to this perturbation 7 . This geometry is 
a compromise, because the magnetic field penetrates our small cluster; however, it lends itself to an extension 
to the full plane, such that only the 4 central plaquettes feel a magnetic field, and the rest of the plane only 



A flux cj> — (f>o can be gauged away, and any physical property, for example the ground state energy, is a 
periodic functions of cf> with period <j)o. 

7 Due to the smallness of the system under study, the hopping parameter tdd that contribute a closed path 
around the central Cu must be kept very small to prevent the magnetic perturbation from growing large and 
destroying the structure of the ground state multiplet. Numerically, the computations were performed with 



D. Superconducting Flux Quantization 
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(24) 
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tdd = ±0.01 eV. 
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FIG. 7.: Four-hole response function R versus the normalized flux, for (a) CU5O4 and (b) CU5O16. The parameter 
values are the same of Table I. The inset shows the geometry of CU5O16 with the closed path. 



experiences a vector potential (see Section III). The relevant response function R is 



and depends on the test flux (j> (in units of (j>o). The results for CU5O4 and CusOie with 4 holes are shown 
in Fig. 7; the same trend is obtained for all clusters. The minimum at 0o/2 is the microscopic precursor of 
the superconducting flux quantization; by looking at A (4) and analyzing the numerical ground state wave 
functions at 4>o/2 one finds that it corresponds to pairs of 1 Ai symmetry which are bound, since the minimum 
at half flux quantum corresponds again to a A(4) < situation. Thus the point symmetry of the wavefunction 
changes from 1 B 2 (x 2 — y 2 ) at <j> = to 1 A 1 (x 2 + y 2 ) at <f> = 4>o/2. 

This finding is a much clearer signature of superconducting flux quantisation than the generally accepted 
presence of the two minima, because it implies that the superconductor remains a superconductor after 
swallowing up the half flux quantum. Therefore, in order to determine the flux dependence of the effective 
interaction, one can compute A(N,<f)). It was found [32] that when A(N,0) < 0, then A(N,cj>o/2) is also 
negative. 

Quite different patterns in which the flux quantization is absent are obtained for two and three holes [29]. 
Also, it is clear that it is IFHubbard that forces the paired holes to screen the vector potential like a charge 2e; 
actually, the non interacting case is quite different [29]. If the input data are modified in such a way that A 
becomes positive by severe distortion of the symmetry, the characteristic central minimum is lost altogether. 

By looking at Fig. 7 we realize that with increasing the flux <j>, the 1 B 2 component of the ground state 
is quickly destabilized by the vector potential, while the 1 Ai component decreases in energy and eventually 
becomes energetically favored in presence of the trapped flux. A similar change of symmetry occurs in 
conventional superconductors where the Cooper wavefunction has s-wave symmetry in absence of magnetic 
flux. 

These findings are required by general symmetry principles. In the absence of tdd, the full invariance Group 
of the cluster is II4 and, as we have seen above, the interacting ground state is degenerate, with 1 A\ and 1 B 2 
components. A nonzero tdd at (j> — reduces the symmetry to the Ci v Subgroup; it turned out numerically 
that with a positive tdd {tdd = 0.01 eV in our calculations) the expectation value of the magnetic perturbation 
is negative on 1 B 2 and positive on 1 A\\ therefore the ground state is 1 B 2 at tdd > but changes symmetry 
if the sign of tdd is reversed. Upon switching the vector potential A, the Cu-Cu hopping is complex and 
chiral, so the symmetry is lowered again from Ci v to its Subgroup Z4, which contains only the rotations. 
Since Z4 is Abelian, there are no degeneracies for a generic <j>, so there are no W = pairs and therefore 
repulsion prevails. With increasing the flux from 0, the ground state energy increases to a maximum. Then it 
decreases because, at (p = (f>o/2, the Cu-Cu hopping of Eq.(24) becomes —tdd, which is real; then the full Ci v 
symmetry is restored, resurrecting the W — pairs. The recovery of d v at <p = (f>o/2 enables us to assign 
the eigenvectors to the irreps, as noted above. The change of symmetry of the pair is also readily understood: 
the perturbation caused by tdd > at <f> — becomes the opposite at half fluxon, so the 1 A\ state is lowest 
now. The signature of superconducting pairing is not only the existence of a well defined second minimum at 
half flux quantum, but also the fact that it corresponds to a pairing (A < ) situation, like at <f> = 0. It is 



R(<t>) 



E (<j>) - E o (0) 



(25) 



tdd 
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gratifying that the present model allows for W — pairs of both symmetries, because this is the only way the 
superconducting flux quantization works. 

Finally, it is worth to note that the Hubbard large U system is stabilized by the intrinsic charge fluctuation. 
Since the mechanism is robust (first order) it is conceivable that these low lying pairing states can survive in 
the vortex core of the underdoped phase [33] . 

III. PAIRING IN THE HUBBARD MODEL: GENERAL THEORY 

The CUO2 plane and the clusters that possess the same Ci v symmetry around a Cu ion have 2-hole eigen- 
states of the kinetic energy with vanishing on-site repulsion (W = pairs). Cluster calculations by exact 
diagonalisation show that these are the quasiparticles that lead to a paired ground state, and have supercon- 
ducting flux-quantisation properties. In this Section, we extend the theory to the full plane, and show that the 
W = quasiparticles are again the natural explanation of superconducting flux quantisation. Moreover, by a 
novel canonical tranformation approach, we calculate the effective interaction W c a between two holes added 
to the ground state of the repulsive Hubbard Model. The method is a particularly efficient way to perform the 
configuration interaction calculation. It is based on the symmetry, and dramatically displays the mechanism 
of pairing in the Cu02 plane. We show that the full configuration interaction calculation can be performed 
recursively. At each step, one decouples a class of virtual excitations while renormalising the matrix elements 
of the kinetic term K and of the Hubbard repulsion ^Hubbard- At the end, one obtains an exact, analytical 
canonical transformation producing an effective Hamiltonian for the dressed pair. In order to get actual num- 
bers, however, we had to neglect the renormalizations in the final formula; this approximation is fully justified 
at weak coupling. Thus, a closed-form analytic expression including the effects of all virtual transitions to 
4-body intermediate states (exchange of an electron-hole pair) is derived. The effective interaction is the result 
of a partial cancellation of positive and negative contributions, so it is not necessarily attractive in all cases; 
the general signature of W = pairs is that the absolute value of the interaction is much smaller than in the 
other cases. Since W = pairs emerge from symmetry alone, they remain W = for any coupling strength 
and are an adequate starting point for a realistic theory. This is the reason why weak coupling expansions 
often provide good approximations at intermediate coupling, as observed by several authors [34] [35] [36] . 

The main scheme is ready to include other interactions which are not considered in the Hubbard Model but 
may be important. We find that the effective interaction is attractive and leads to a Cooper-like instability of 
the Fermi liquid, while it is repulsive for triplet pairs. From W e s, an integral equation for the pair eigenfunction 
can be derived; the binding energy |A| of the pairs is in the range of tens of meV, which is not comparable 
to any of the U and t input parameters. The reason is that the interaction, which vanishes identically for the 
bare W = pairs, remains dynamically small for the dressed quasiparticles. This suggests that a weak coupling 
theory may be useful to study the pairing force, despite the fact that U is not small compared to t. 

A. Band Structure 

Let us consider the three-band Hubbard Hamiltonian. The kinetic term is 

K = ^ t ss ,cl a c s , a + ^2 £sc\ a c sa (26) 

(s,s'} s 

where t 3S i = t pd if s and s' represent a Copper site and a nearest neighbor Oxygen site, or viceversa, and zero 
otherwise, while e s — e p if s is an Oxygen site and zero otherwise. It is convenient to identify each unit cell 
by its cartesian coordinates r = (i x ,i y ), i x ,i y = 1, . . . , Na (N% = |A|/3 being the total number of unit cell) 
in such a way that r, r + a and r + b are the positions of the Cu and of the two Oxygens, respectively. The 
one-body eigenstates of K are 

Kc^f\0)=e^c^\0), k=(k x ,k y ) = ^-(i x ,i y ), i x ,i y = l,...,N A , (27) 

where v is a band index and periodic boundary conditions are chosen. According to Bloch's theorem they can 
be written as 

=5>LV> M (k, S ), ^(k, S )=e lk V ) (k, S ) (28) 
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where tp("' has the lattice periodicity and r is the position of the cell in which the site s lies. The non-bonding 
band is characterized by ip ( - nb ' 1 (k, 0) = 0; one then finds 



^" 6) (k,a) 



cos [^] 



iV A 



(cos 2 [^] + cos 2 



¥><" k) (k,b) 



N A 



(COS 2 [!f] + COS 2 ) 



For the other bands, one obtains the eigenvalue equation 



(± ) _ e p ±r 



pd 



+ cos 2 



(+ for the antibonding and - for the bonding band) and the Cu amplitudes 



^ (±) (k,0) = -^- 



e(±) 



iV A V (2e (±) - e p ) ' 
while the O amplitudes are obtained from 

( £ p _e k ±) ) V (±) (k,a) +2t pd cc 



^ (±) (k,0)=0, 



and 



(e P - e^) V (±) (k,b) + 2t pd cos 



y> (±) (k,0) = 0. 



To write VKijubbard in terms of the c^J we must invert Eq.(28): cJ CT = ^ kl/ (k, s). Thus, 

ki . .k4 v\..i/^ 

Using a shorthand notation, where the band indices are understood, 

U (kl, k 3 , k 2 , k 4 ) = ^ e -(kl+k3-k 2 -k 4 ). r ^ Utip . s) ^ (k3) s) ^ (k2j s) ^ fi) ^ 

r sfEccll 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



where X/seccii mv °l ves a sum over a single cell. The umklapp processes involving the basis vectors of the 
reciprocal lattice produce a minus sign on one of the Oxygens. The matrix elements (35) are proportional to 
iVX 2 , like any bona fide two-body interaction. 

The one-band Hubbard Model is even more elementary; we postpone the details to the next Section. 



B. W — pairs and Superconducting Flux Quantization 

In this Section we demonstrate that W = pairs with vanishing momentum and formed by holes in 
degenerate states exist whatever is their energy. Thus, in the thermodinamic limit, there are no constraints 
on the occupation number; any even number is "magic". Instead of using the W — theorem and the full 
symmetry, that would in principle yield all the W = pairs, we follow the simpler route of projecting a single 
determinantal state on the irreps of Gav [37] [38] . This partial use of the symmetry still gives enough solutions 
to demonstrate pairing. 

Let us first consider the Hubbard Model with three bands. Omitting band indices again, we shall mean 

|d[k])=4 iT cl kil |0> (36) 

to be a two- hole determinantal state derived from the k eigenfunctions. Since tp(— k, s) = <p (k, s), which is 
required by time-reversal symmetry, the combination |d[k]) + \d[— k]} is singlet and |rf[k]} — \d[— k]) is triplet. 
We introduce the determinants |d[i?k]) = |d[kfl]), R G Cav, and the projected states 

|*„[k]> = -^ Y, X W (fi)l« (37) 

Rec 4v 
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FIG. 8.: (A) Pattern of the vector potential A due to 4 flux tubes (empty dots) carrying flux <f>. X stands for Cu. 
The line integral of A along each bond parallel to the arrow is 0/2. (B) The geometry of the torus and of the magnetic 
field used to discuss the flux quantisation properties of the W = pairs. 



where x (R) is the character of the operation R in the irreducible representation rj. In the non-degenerate 
irreps, the operations that produce opposite ki? have the same character, and the corresponding projections 
lead to singlets. Let Ri, i — 1, ..8 denote the operations of Cav and k, k' any two vectors in the Brillouin Zone 
(BZ). Consider any two-body operator O, which is symmetric (R\ORi = 6), and the matrix with elements 
Oi,j = (d\k]\R\ORj ;|d[k']), where k and k' may be taken to be in the same or in different bands. This matrix 
can be diagonalized by Group Theory alone. Indeed, for each irrep rj, consider the normalized vector with 
components -^gx^ (Ri)', this is an eigenvector of the operator matrix R\ORj, as we can check by noting that 

f Ei j X M (Ri) R]ORjX^ (Rj) = OP M2 S(rj, ??'); hence, RjORj is diagonal on the basis of the symmetry 

projected states. The square of the projection operator is P^ — V&P^. Now, taking the matrix element 
between the k and k' determinants, we get the eigenvalues in terms of determinantal matrix elements: 

0(7j,k,k') =]T X (,7) (.R)OK(k,k') (38) 

R 

where Or (k, k') = (d[k|0|d[k^.]). If O is identified with Wriubbard, the determinantal matrix element does 
not depend on the sign of the components of k', because the <p's depend on k only through a cosine. The only 
thing that matters is that , o+ and a- exchange k' x and k' y , while the other operations leave them 

in place. Consequently, (W Hu bbard)i = (W / Hubba r d)c 2 = (W / Hubbard)<r x = (WHubbard)^ and (W / Hu bbard) r ,(+) = 

(WHubbard) c (-) = (W / Hubbard)<r + = (^Hubbard) a_ ■ Therefore 

WWbbard (M 2 ) = WW-bard (^2) = 0. (39) 

These are W — pairs, like those studied previously [29], but with an important change, since in the cluster 
calculations the symmetries of W = pairs were found [26] to be 1 £>2 and Mi. The reason for this change is 
a twofold size effect. On one hand, 1 Ai pairs have the W — property only in the small clusters, having the 
topology of a cross, and belonging to the LI4 Group, but do not generalize as such to the full plane, when the 
symmetry is lowered to C'4 V ; on the other hand, the small clusters admit no solutions of 1 A2 symmetry at all 
if only degenerate states are used. One necessary condition for pairing in clusters is that the least bound holes 
form such a pair, and this dictates conditions on the occupation number. In the full plane, however, W = 
pairs exist at the Fermi level for any filling. 

Let us now explore what happen if a magnetic field is turned on [32]. Consider the pattern of Figure 8 (A). 
Here, the Copper sites are marked by X and the Oxygen sites by O; the black dots stand for tubes carrying flux 
4> each, symmetrically disposed around the central Cu. Varying <j) by an integer multiple of <f>o corresponds to 
a gauge transformation leaving all the physical properties invariant. The arrows help to visualise a convenient 
choice of the vector potential at general <f>. Namely, running along an oriented bond in the sense of the arrow, 

f A-dr=t (40) 
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along the other Cu-O bonds, not marked in the figure, / A-dr = 0. One sees that in this way the flux through 
any closed path corresponds to the number of tubes surrounded by the path. The reflection operations of 
dv are equivalent to <j> — > — 0, reverse the directions of the arrows and for a generic <j> the symmetry Group 
reduces to Z4. However, at 4> = 4>o/2 the reversal of the magnetic field in the tubes corresponds to a jump 
by <f>o, and this is equivalent to a gauge transformation: this implies that the symmetry Group gets bigger, 
the new symmetry operations being reflections supplemented by a gauge transformation. Indeed, it follows 
from Eq.(40) that the hopping parameter becomes it p d along the arrows, while it remains equal to t p d along 
the unmarked bonds of Figure 8 (A). Any reflection operation simply changes the signs of all the hoppings 
along the marked bonds. Now consider the unitary transformation U which changes the signs of all the Cu 
orbitals along both diagonal, except the central Cu. Since U also has the effect of reversing all the arrows, 
a ■ U is a symmetry, for all reflections a in Ca v . Moreover, since the product of two reflections is a rotation, 
the Group Ca v including the rotations and the reflections multiplied by U is isomorphic to C\ v . The W = 
pairs appropriate for half a flux quantum must involve two holes belonging to the degenerate irrep of Ca v . In 
this way, at 4> = 4>o/2 the full symmetry is restored, allowing again for pairing and negative A. The W = 
quasiparticles have just the correct symmetry properties in the presence of the vector potential to provide 
superconducting flux quantisation in macroscopic systems. 

W = pairs exist only in planar lattice structure with a certain point symmetry Group. So we expect that 
the CUO2 geometry is not the only one where W = property is present. In the following we want to show 
that in the one-band Hubbard Model too, where the point symmetry Group is still Cav , we can build W = 
pairs and that they continue to exist even in the presence of magnetic field if and only if the magnetic flux is 
quantized in units of <j>o/2 [37]. 

Let us start with the one-body Hamiltonian defined on a square lattice with N\ x N\ sites 8 (for the moment 
we don't care of the spin index), 

K = \ ^ (4cr' + h-c) + e cjc r , r = (i x ,i v ), i x ,i v = 1, . . . , N A (41) 

(r,r') r 

where (r, r') stands for nearest neighbors. We take the x,y axes parallel to the bonds and fold the plane 
on itself after N\ lattice spacings in both directions, forming the torus shown in Figure 8 (B). Since we are 
dealing with complex fermion fields, the most general boundary conditions are 

c r+Ta = e a c T c r+Tb = e^c v (42) 

where r a and r b are the torus periods: r a — (Na,0) and Tt = (0, Na). The fermion fields can then be expanded 
as 

k 

with t = (o/Na, /3/Na)- Turning on a costant magnetic field B = (B a , Bb) parallel to the surface of the torus, 
the flux through a circular surface S a (E 6 ) of radius ro = Na/2-k perpendicular to the 6 (a) direction will be 

<I>t = B nda=(h A ■ dl = 2nr A T , r = a,b (44) 

with n the surface versor. The above equation enables us to compute, using the Peierls prescription, how the 
presence of this particular constant magnetic field changes the original Hamiltonian K of Eq.(41). The final 
result is 



if — > - 

2 



(e^^ /M clc r ,+h.c.)+ (e^ /M clc r ,+h.c.) 

(r,r'}a <r,r'}i, 



c\c v (4£ 



where (r, r') T are two nearest neighbor sites along the r = a,b direction. The spectrum of the new Hamiltonian 
is the same of the original one if a — > a + 2n(f> a /(f>o and f3 — > ft + 2n(f>b/4>o- In this sense we can say that a 



8 Although we know that the actual cuprate superconductors present a band structure with three dominant 
bands, the one-band Hubbard Model tries to mimic the presence of the charge transfer gap A p d by means of 
an effective value of the Coulomb repulsion and thus it presents only two bands. The "Oxygen" band becomes 
the lower Hubbard band of this model. 
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change in boundary conditions is the same of turning on gauge fields [39] . Of course the choice which physically 
corresponds to the unperturbed infinite plane is a — ft = (Born- von Karman boundary conditions). 
Denoting with 

e ia 

□ 4> b (46) 

<j)a 

the configuration for which the boundary conditions are e la and e l13 and the flux is 4> a and (pi along the the 
a and b direction respectively, explicit calculations show that for the Hamiltonian (45), W = pairs builded 
from determinantal states of the form in Eq.(36) of 1 A2, 1 £?2 and 1 Bi symmetry exist if and only if one of the 
following configurations is realized 

e ia 

e l/3 □ (np - (47) 

(n a - a/n)<p /2 

with n a and np integers such that (— l) n ° = (— l)™* 3 . 

From Eq.(47) we see that only two configurations allow for W — pairs in absence of magnetic field; they 
correspond to choose boundary conditions both periodic or antiperiodic along two orthogonal directions. With 
such boundary conditions W = pairs exist even in a magnetic field if and only if the flux is quantised in 
multiples of half the fundamental fluxon. For any other choice of boundary conditions, flux is still quantised 
in a similar manner, but we need a zero point flux in order to have W = pairs, so that they do not exist 
at zero magnetic field. If we assume dressed W — pairs to be quasiparticles of the superconducting phase, 
then flux is quantized in a superconducting way. In the next Section we will prove rigorously that this is the 
case, by studying the role of W = pairs in the many-body problem and showing how they become a genuine 
bound state by means of the electronic correlation. 



C. Canonical Transformation 



In the Cooper theory, an effective interaction involving phonons is introduced via an approximate canonical 
transformation. In the Hubbard Model, the holes can exchange particle-hole pairs, and more complex electronic 
excitations, rather than phonons, but one can approach the problem in a similar way [37] [38] [40]. 

Suppose the system is in its ground state with Fermi energy ef and a couple of extra holes are added. 
We wish to show that by a canonical transformation one obtains an effective Hamiltonian which describes 
the propagation of a pair of dressed holes, and includes all many-body effects. Including all the diagrams 
of a generalized RPA would lead to something like the FLEX approximation [41] whose implications for the 
superconductivity in the three-band Hubbard Model have been explored recently in a series of papers [42]. 
A related self-consistent and conserving T-matrix approximation has been proposed by Dahm and Tewordt 
[43] for the excitation spectra in the 2d Hubbard Model; we mention incidentally that recently diagrammatic 
methods have been successfully applied to the photoelectron spectra of the cuprates in other contexts too, like 
the spin-fermion model [44]. 

The exact many-body ground state with two added holes may be expanded in terms of excitations over the 
vacuum (the non-interacting Fermi sphere) by a configuration interaction: 

|*o(JV + 2)) = ^a m |m) + ^a a |a) + ^a fj \p) + .... (48) 

maf3 

here m runs over pair states, a over 4-body states (2 holes and 1 e-h pair), (3 over 6-body ones (2 holes and 
2 e-h pairs), and so on 9 . To set up the Schrodinger equation, we consider the effects of the operators on the 
terms of |* (JV + 2)). We write: 

K\m)=e m \m), K\a)=e a \a), K\0) = e p \0), ... (49) 



9 Many configurations contribute to the interacting ground state and we need a complete set to expand it 
exactly. However, as long as the set is complete, the expansion is exact and well defined only if the non- 
interacting ground state is non-degenerate. We assume this is the case in Eq.(48). Besides, we stress that 
Eq.(48) is configuration interaction, not a perturbative expansion. 
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and since W^Hubbard can create or destroy up to 2 e-h pairs, 



(50) 



W m i^ m does not contribute to the effective interaction for W = pairs in our model; however we keep it for 
the sake of generality. It allows to introduce the effect of the exchange of phonons and other quasiparticles 
that we arc not considering. For clarity let us first write the equations that include explicitly up to 6-body 
states; then we have 



WWbbardla) = ^V m , a |m) + ^V a /,«|a') + 



(51) 



where scattering between 4-body states is allowed by the second term, and 

^Hubbard!/?) = Y, W ™'>P K) + E^ ^ + E^ \P) 



(52) 



In principle, the Wpi t p term can be eliminated by taking linear combinations of the complete set of f3 states 
such that (if + lVHubbard)/3',^ = e'pSppr. when this is done, we get a "self-energy" correction to ep and a mixing 
of the vertices, without altering the structure of the equations. The Schrodinger equation yields equations for 
the coefficients a m , a a and dp 



(e m - B ) a m + ^W m , m /a m / + ^W m , a o a + ^W m ,pap = 



(53) 



(e a - -Bo) a a + ^W a , m /a m / + ^W aiQ 'a a / + ^Wq^pap = 



(54) 



(e'p - E ) ap + y^Wp, m >a m > + /~~]Wp, a 'a a ' = 



(55) 



where Bo is the interacting ground state energy. Then, we exactly decouple the 6-body states by solving the 
equation for ap and substituting into Eqs.(53)(54), getting: 



(e m — i?o) »m + 



W 



m,m' 



+ 



E 



W mt pWp >m > 
E - e'„ 



W m , a + ^ 

13 



W m ,pWp„ 
E - £% 



a a =0 (56) 



(e Q - Bo) a a + 



w , + V 

Y¥ a.m' \ 7 

13 



W a ,pWp, 
Bn — e' a 



V + 



E 



-r^W^pWp, 
^ E -e' 



a„i = 



(57) 



Thus we see that the role of 6-body states is just to renormalize the interaction between 2-body and 4-body 
ones, and for the rest they may be forgotten. If Bo is outside the continuum of excitations, as we shall show 
in the next Section, the corrections are finite, and experience with clusters suggests that they are small. Had 
we included 8-body excitations, we could have eliminated them by solving the system for their coefficients 
and substituting, thus reducing to the above problem with further renormalizations. In principle, the method 
applies to all the higher order interactions, and we can recast our problem as if only 2 and 4-body states 
existed. Again, the W a i >a term can be eliminated by taking linear combinations of the a states: when this 
is done, we get a "self-energy" correction to e a and a mixing of the W m ,a vertices. Denoting W' and e the 
renormalised quantities, which can be determined by the procedure outlined above, one can solve for a a . The 
eigenvalue problem becomes 



(Bo- a„ 



(58) 



where 



(m\S[E ]\m')=J2 



[m\W'\a){a\W'\m'} 
Bo — e' a 



(59) 
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Eq.(58) is of the form of a Schrodinger equation with eigenvalue Eo for pairs with an effective interaction 
W' + S. Then we interpret a m as the wave function of the dressed pair, which is acted upon by an effective 
Hamiltonian impair. This way of looking at Eq.(58) is perfectly consistent, despite the presence of the many- 
body eigenvalue Eo, because we are not compelled to reference the energy eigenvalues to the hole vacuum. We 
note that if we shift the kinetic operator K by an arbitrary constant AE in Eq.(49), by setting K' = K — AE, 
the same shift applies to the eigenvalues e m ,e a , tp and so on, and also to the renormalised quantities like e' a , e'p. 
Therefore, the matrix elements of the scattering operator S are unaffected by the shift. Thus we can reference 
Eo to a new energy origin by shifting the diagonal terms in Eq.(58) without changing the off-diagonal terms. 
Since we wish to regard Eq.(58) as a Cooper- like equation for the pair, it is natural to set AE equal to the 
interacting ground state energy eigenvalue for the N hole system, relative to the hole vacuum. The energy of 
two independent holes, relative to the N background, is 2e^; when the effective interaction is accounted for, 
the energy of the two bound holes is 2ef + A, where | A| is the binding energy. 

The change from the full many-body i/Hubbard to impair is the canonical transformation to dressed pairs 
we were looking for. In the new picture, the holes interact through an effective vertex with infinitely many 
contributions. The linked contributions represent repeated exchange of electron-hole pairs, and may contain 
self-energy insertions; all these contributions make up the effective interaction. The unlinked diagrams are 
pure self-energy. Thus, the scattering operator S is of the form S = W e g + F, where W e s is the effective 
interaction between dressed holes, while F is a forward scattering operator, diagonal in the pair indices m ,m' 
which accounts for the self-energy corrections of the one-body propagators: it is evident from Eq.(58) that it 
just redefines the dispersion law e m and we suppose it be included in the kinetic term K. 

So, letting \a) — ~^2 m a m \m), the effective Schrodinger equation for the pair reads 

(K + W' + W eS )\a) = E \a) (60) 

and pairing takes place if Eo = 2ef — |A|; then |A| is the binding energy of the pair. Any other pairing 
mechanism not considered here, like off-site interactions, inter-planar coupling and phonons, can be included 
as an extra contribution to W' m i m which just adds to W e g. 

We emphasized the fact that in principle the canonical transformation is exact because in this way the 
framework does not require U jt to be small. However to test the instability of the Fermi liquid towards 
pairing it is sufficient to study the amplitudes a m of the m states. In the weak coupling limit this can be done 
truncating the expansion (48) to the a states because, as we have shown, the inclusion of the f3, 7, . . . states 
produces a renormalization of the matrix elements of higher order in U, leaving the structure of the equations 
unaltered 10 . 

Let |$o(A r )) be the non-interacting ground state formed by filling all the energy levels up to €f(N) for both 
spin up and down. If the added pair has vanishing total momentum, like the one in Eq.(36), the a states are 

I") = 4< +q+ k 2 ,T c k2,|C t _ qa c! k , J$oW) (61) 

where Ck 2 ,| creates an electron state; hence £k 2 < £f < £k'+ q +k 2 , e i, e k'- The a states with opposite spin 
indices arise from the specular diagram, contribute similarly and yield a factor of 2 at the end. The interaction 
matrix element is 

(<I>o(iV)|c_ k / >i c_ q>i 4 2a c k / +q+k2iT W 7 Hubba r d|d[k]) = 
<S(q - k)C/(k' + q + k 2 , -k', k, ka) - 8{k' - k)[/(k' + q + k 2 , q, k, k 2 ) (62) 

where the U matrix is computed on the Bloch basis, see Section III A. This expression holds for any lattice 
with translational invariance. Working out Eq.(59) we find that the product in the numerator yields 4 terms; 
two are proportional to <5(p — q) and belong to F, while the cross terms yield identical contributions to W e ff- 
Thus we obtain the following effective interaction between the members of W = pairs: 



10 We want to stress that the truncated expansion of |^o (iV + 2)) in Eq.(48) doesn't give a good approximation 
of the interacting ground state wave function but only of its weak coupling a m amplitudes. Similarly in the 
BCS Model, from the Cooper equation (obtained by truncating |^o(^V + 2)) to the m states) we can estimate 
only the pair coefficients of the ground state and not its full structure. Nevertheless this is enough to study 
bound states formation; indeed the energy gap of the pair in the Cooper theory and in the many-body BCS 
theory are identical. Moreover, the pair wave function of the Cooper Model resembles the one of the many-body 
BCS Model [45]. 
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N A 


p = N/N'i 


—A (meV) 


V cH (eV) 


-A asympt (meV) 


18 


1.13 


121.9 


7.8 


41.6 


20 


1.16 


42.2 


5.0 


9.0 


24 


1.14 


59.7 


7.0 


28.9 


30 


1.14 


56. 


5.7 


13.2 


40 


1.16 


30.5 


6.6 


23.4 



TABLE VI.: Binding energy of 1 £?2 pairs in supercells. 
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while it can be shown that the corresponding matrix elements of the forward scattering operator F vanish 
identically as in the fully symmetric clusters [46]. The sum in Eq.(63) is over occupied k with empty p + qa + k. 
Note that W e g does not depend on the sign of U. The diagonal elements of Eq.(63) are also clearly related 
to the A (2) expression (20) derived from perturbation theory for the fully symmetric clusters [29], which are 
special cases of the present theory. The expressions (20) and (63) are characterized by the symmetry-induced 
quantum mechanical interference of several terms. The sum can be positive or negative depending on the 
irrep. This interference produces a partial cancellation, and the absolute value of the result is typically much 
smaller than individual contributions. This explains why the interaction is dynamically small for W = pairs: 
they have no direct interaction, and because of the interference the effective interaction is reduced compared 
to what one could expect by a rough order-of-magnitude estimate. However, the presence of the 6 functions 
and the anisotropy of the integrands prevent a total cancellation. 

In Eq.(20), the vanishing of the denominators is prevented by the discrete spectrum of the cluster. In the 
thermodynamic limit this complication disappears in the full expression (63) since we are interested in bound 
pairs; then, Eo is below the continuum and the denominators never vanish. 

The p and q indices run over 1/8 of the empty BZ for W = pairs. We denote such a set of empty states 
-D(£f)/8, and cast the result in the form of a (Cooper-like) Schrodinger equation 

2e k a k + Y W * s ( k ' k a k' = £ offlk (64) 

k'£D(e F )/8 

for a self-consistent calculation of Eo (since W e ff depends on the solution). 

D. Pairing in the CuC>2 Plane: Numerical Results 

Using the analytical expression for the effective interaction in the full plane, Eq.(63), we have performed 
exploratory numerical estimates of A by working on supercells of Na x Na = |A|/3 cells, with periodic 
boundary conditions. For the sake of simplicity, we shall neglect the minor contributions from the higher 
bands and consider the dominant intra-band processes, in which empty states belong to the bonding band. 
Here we solved the problem in a virtually exact way for Na up to 40. Several good supercell calculations 
devoted to the problem of pairing have been reported to date [7], but no conclusive evidence was reached, 
because of the difficulty of dealing with size effects. First, we searched for triplet solutions with negative 
energy without success, since, as in the clusters, W e s is repulsive for triplets. On the contrary, W = singlets 
did show pairing, in line with our previous findings in small clusters [18] [19] [23] [29] . We took as input data 
the set of current parameters, already used for clusters, that is (in eV) t p a = 1.3, e p — 3.5, e<j = 0, U p = 6s, 
Ud = 5.3s, where s is a scale factor induced by renormalization. Since screening excitations are explicitly 
accounted for in the Hamiltonian, it is likely that U is a bare (unscreened) quantity, which justifies s > 1. 
A stronger interaction causes smaller pairs and speeds up convergence within attainable supercell sizes. In 
Table VI, we report the results for 1 i?2 pairs at s = -J= = 2.121 with ef = —1.3 eV (half filling corresponds 
to £f = —1.384 eV). We see that although |A| decreases monotonically with increasing supercell size, V e s, see 
below for its definition, is not dropping to zero, but remains fairly stable around 6 to 7 eV. With supercell 
sizes Na > 40 calculations become hard. Since we are concerned with the asymptotic behavior for Na — > oo 
and A depends on U's and Na in a peculiar way, increasing with Na for large U's and dropping for small 
IPs, we need a simple solvable model in supercells and in the infinite plane to make reliable extrapolations of 
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e F (eV) 


lo d (eV) 


p F (eVVcell) 


-1.35 


0.4545 


0.0486 


-1.3 


0.492 


0.0415 


-1.2 


0.520 


0.0338 


-1.1 


0.5237 


0.02905 



TABLE VII.: Effective Cooper parameters for fitting the V e g dependence of A at several ep values. 
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6.9 


135. 
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27. 


6.3 


198. 


13.8 


460. 20.8 


00 


8.8 


6.3 


131. 


13.8 


333. 20.8 



TABLE VIII.: Data for 1 B 2 pairs at e F = -1.2 eV; V cS is in eV. 



numerical results. To this end, we define the Uniform Interaction Model (UIM) in which a constant negative 
interaction — Vff, V e g > prevails for k and k' in Z?(ef)/8. In the Cooper formula for the binding energy of 
the pair one has 

|A| = (65) 
e ~i PF —i 

where pf is the density of states at the Fermi level and —7 is a (constant) attractive interaction per cell in 
a shell of thickness 2lod surrounding the Fermi surface (the orders of magnitude are lod ~ 25 meV, pF ~ 0.3 
eV _1 /cell and 7 w 1 eV/cell, and | A| is in the meV range); on the other hand, in the UIM no Debye frequency 
ujd is involved, and all the empty states in D(tF)/8 contribute to A. Both the Cooper and the UIM models 
can be solved in supercell calculations, and show the same qualitative behavior with increasing Na- 

For small Na, A is of the same order as V e ff, while in the thermodynamic limit (Na — > 00), the asymptotic 
value of A, A asymp t, can easily be estimated. Setting W e s = —Vbs/Na, we get from Eq.(64) 

(2e k - E )a k = J^y j dk ' a k' = const. (66) 

Writing a k = const/(2e k — Eo) and Eo = 2cf — |A|, we obtain 

8 f° dep (e) 



V e g J 2 (e — Ef) + I A asymp t I 



(67) 



Here, p is not a constant, as in the Cooper formula, and its values in the integration range are sensitive to 
filling; howewer its order of magnitude does not differ much from the BCS case except at half filling. Good 
fits can be obtain using the Cooper formula with effective lod and pf parameters, like those shown in Table 
VII. 

In Table VIII, we consider 1 B2 pairs at = —1.2 eV and report values of A and V e ff for supercell 
calculations at Na=12, 20 and 30 for 3 values of the scale factor s. The filling is p ~ 1.3 at Na = 30. Fixing 
Vff at the value of Na = 30 we calculate A asymp t for Na — > 00. 

We see that the relatively mild Na dependence of V e g supports the use of the UIM to extrapolate the results 
to the thermodynamic limit, and there is a clear indication of pairing with sizable binding energies. The s 
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-A (meV) Vefi 


12 


163. 


9. 


409. 
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20 
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371. 


17.6 


635.5 24.3 


30 


80. 


10.3 


270. 


17.66 


520. 24.9 


00 


35.6 


10.3 


170.8 


17.66 


362. 24.9 



TABLE IX.: Data for 1 i?2 pairs at ep = —1.1 eV; V ctt - is in eV. 
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dependence of V e g is roughly linear, while A depends exponentially on s. Table IX presents the results for 
1 f?2 pairs at ef = — 1.1 eV (p w 1.4 at Aa = 30). The trend is similar, but V e g is seen to increase with doping. 

The results for the 1 A2 pairs are seen to lead to bound states as well, with comparable A values [37]; the 
trend with doping is opposite, however, and the gap is nearly closing at ef = — 1.1 eV. A necessary condition for 
superconducting flux quantization is that two kinds of pairs of similar binding energy and different symmetries 
exist [26]. A similar conclusion was reached independently by other authors [47]. Moreover, evidence of mixed 
(s + id) symmetry for the pairing state has been amply reported in angle- resolved photoemission studies [48]. 
This remark leads to the prediction that in this model superconductive pairing disappears with increasing p, 
while a different sort of pairing prevails; in reality the CuC>2 plane prefers to distort at excessive doping, and 
in a distorted plane the present mechanism, based on symmetry, could be destroyed. 



E. Antiferromagnetic Ground State at Half Filling 

Up to now, we learned how to master the situation when the number N of holes in the A/a x A?a system is 
such that the non-interacting ground state |<J>o(AT)) is a single non-degenerate determinant (the Fermi surface 
is totally filled). Then, switching the repulsion Wmibbard on, we know how to calculate the effective interaction 
W e s between two extra holes added to the system by the expansion (48) . However we are particularly interested 
in the doped antiferromagnet, and the antiferromagnetic ground state occurs at half filling, not in a filled-shell 
situation. 

We want to study the doped antiferromagnet since there are strong indications that the Fermi liquid is 
unstable towards pairing near half filling; they come from diagrammatic approaches [41], renormalization 
group techniques [49] [50] and also cluster diagonalizations [29] [51] [52]. Therefore exact results on the half 
filled Hubbard Model may be relevant to antiferromagnetism and to the mechanism of the superconducting 
instability as well. On the other hand, pairing by the "W = 0" mechanism does not require that the pair be 
created on a single-determinant state. 

For simplicity in the rest of this paper we concentrate on the one-band Hubbard Model, defined by Eq.(41) 
plus the usual on-site Hubbard interaction WHubbard: 

^Hubbard = K + V^Hubbard = | ^ (4 CT C r ' CT + h.C.) + U ^ n r} n ll (68) 

{r,r')cr r 

with r = (i x ,i v ), i x ,i y = 1, . . . , N\, the sum on (r, r') is over the pairs of nearest neighbors sites and h va 
is the number operator on the site r of spin a 11 . The point symmetry is Ci v \ besides, /^Hubbard is invariant 
under the commutative Group of translations T and hence the space Group [53] G — T ® C'4 V ; ® means 
the semidirect product. In terms of the Fourier expanded fermion operators (periodic boundary conditions) 
Ck = ^2 r e lk r c r , we have K = ^ k ekcj^Cko- with e k = 2t[cos k x + cos k y ]. Then, the one-body plane wave 
state cj^lO) = jka) is an eigenstate of K. 

The 4x4 cluster is half filled at A^ = 16. Below, we shall consider its ground state with N = 4 [54] and 
N = 14 [51] [55] holes. For U — > 0, A" = 4 means that 2 holes are added to a closed shell, while N = 14 
corresponds to 2 electrons added to the half filled system. We demonstrate that in both cases the structure of 
the ground state must be interpreted in terms of pairing 12 . However the non-interacting ground state at half 
filling is degenerate and in order to deal with N — 14 we must extend the above formalism. 

As a preliminary measure in this Section we build the exact ground state of the Hubbard Model at half 
filling and weak coupling for a general even Na [56] [57]. Once a unique non-interacting ground state is 
determined, one can use the non-perturbative canonical transformation to test the instability of the system 
towards pairing; this will be done in the next Section for Na = 4. 

The starting point is the following property of the number operator n r = cjc r (for the moment we omit the 
spin index). 

Theorem : Let S be an arbitrary set of plane-wave eigenstates { | ) } of K and {n r )ij = (ki|n r jkj) = 
1 e »(k»-k 3 )-r j.jje matrix of h T in S. This matrix has eigenvalues Ai = jjSj- and A2 = . . . = \\s\ = 0. 

A A 



11 Here we have dropped the one-body on-site potential e^ro- c raC r <T, since we work with a fixed number of 
spin up and down particles. 
12 Hole pairing in the former case and electron pairing in the latter. 
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Note that |<S| < N/[; if \S\ = N\ the set is complete, like the set of all orbitals, and the theorem is trivial (a 
particle on site r is the n r eigenstate with eigenvalue 1). Otherwise, if \S\ < N%, the theorem is an immediate 
consequence of the fact [57] that 

det|(n r ) l3 -A5 l3 | = (-A) |s| - 1 (i^-A), Vr. (69) 

A 

Let Shf denote the set (or shell) of the k wave vectors such that — 0. At half filling (A^ particles) for 
U — the Shf shell is half occupied, while all k orbitals such that £k < are filled. The k vectors of Shf 
lie on the square having vertices (±7r,0) and (0, ±7r); one readily realizes that the dimension of the set Shf is 
\Shf\ = 2Na — 2. Since Aa is even and .^Hubbard commutes with the total spin operators, at half filling every 

/ 27V 2 \ ^ 

ground state of K is represented in the S z = subspace. Thus, K has [ jy ^ j degenerate unperturbed 

ground state configurations with S z = 0. We wish to study below how this degeneracy is removed by the 
Hubbard interaction WHubbard already in first-order perturbation theory. Actually most of the degeneracy 
is removed in first-order, and with the help of a theorem by Lieb [58] we shall be able to single out the 
true, unique ground state of ^Hubbard in Eq.(68) 13 . The first-order splitting of the degeneracy is obtained 
by diagonalizing the WHubbard matrix over the unperturbed basis; like in elementary atomic physics, the 
filled shells just produce a constant shift of all the eigenvalues and for the rest may be ignored in first-order 
calculations. In other terms we consider the truncated Hilbert space H spanned by the states of Aa — 1 holes 
of each spin in Shf, and we want the exact ground state(s) of WHubbard in W; by construction TL is in the 
kernel of K , so the ground state of WHubbard is the ground state of -^Hubbard as well. We call W = state any 
vector in 7i which also belongs to the kernel of WHubbard. Since the lowest eigenvalue of WHubbard is zero, it 
is evident that any W — state is a ground state of -ff H ubbard- We want to calculate the unique ground state 
of the Hubbard Hamiltonian for U = + at half filling which is a W = singlet with 2Na — 2 holes in Sh / 
(filled shells are understood) 14 . 

To diagonalize the local operator WHubbard in closed form we need to set up a local basis set of one-body 
states. If Shf were the complete set of plane- wave states k, the new basis would be trivially obtained by a 
Fourier transformation, but this is not the case. The question is: how can we define for each site r the local 
counterparts of k states using only those that belong to a degenerate level? The answer is: build a set {|^L r) » 
of orbitals such that the number operator n r and the Dirac characters of the point symmetry Group Cav are 
diagonal. Using such a basis set for the half-filled shell the unique properties of the antiferromagnetic ground 
state become simple and transparent. The eigenvectors \<p^} of n T =o and those \y>^) of other sites r are 
connected by translation and also by a unitary transformation, or change of basis set. Picking r = ei, I = x 
means e; = (1, 0) or transfer by one step towards the right and I = y means ei = (0, 1) or transfer by one step 
upwards. The unitary transformation reads: 

2JV A -2 2JV A -2 

i^ !) >= E wT)(*T\^ i) )= E i^ 0) >tw (7i) 

0=1 13=1 

The transfer matrix Ti knows all the translational and point symmetry of the system, and will turn out to be 
very special. 

For large Aa, to find {Iv 3 ^}} it is convenient to separate the k's of Shf in irreducible representations of the 
space Group G= Ca v ® T. Choosing an arbitrary k £ Shf with k x > k y > 0, the set of vectors 7?;k, where 
Ri £ Ci v , is a (translationally invariant) basis for an irrep of G. The high symmetry vectors (0, 7r) and (tt,0) 
always trasform among themselves and are the basis of the only two-dimensional irrep of G, which exists for 
any Aa- If Aa/2 is even, one also finds the high symmetry wavevectors k = (±7r/2, ±7t/2) which mix among 



The Lieb theorem states that at half filling the ground state for a bipartite lattice is unique and has spin 
||<5>i| — |<S2||/2 where |<Si| ( | <S*2 1 ) is the number of sites in the Si (S2) sublattice. 
14 Equivalently, we shall find the exact (2Aa — 2)-body ground state of the effective Hamiltonian 

tfeff=-^4 E 5 ( k l+ k 2- k 3-k4)4 lT C k 2 | C k3l C k 4 T (70) 

A k 1 ,k 2 ,k 3 ,k 4 eS h/ 

where 5(G) = 1 if G is a reciprocal lattice vector and zero otherwise. 
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themselves and yield a four- dimensional irrep. In general, the vectors J?ik are all different, so all the other 
irreps of G have dimension 8, the number of operations of the point Group Cav 

Next, we show how to build our local basis set and derive W = states for each kind of irreps of G. For 
illustration, we will consider the case Na = 4; then Shf contains the bases of two irreps of G, of dimensions 
2 and 4. The one with basis = (tv, 0), k_g = (0, n) breaks into A\ © B\ in Cn v . 

The eigenstates of (n r =o)ij = (kijn r= o|kj), with i,j — A, B, are \tp' Al ) = -^(|kA) + |ks)) with Ai = 1/8 and 
= "^Tf (I^a) — |ks)) with A2 = 0. Since under translation by a lattice step Ti along the I — x,y direction 

|k) -> e lk > |k), using Eq.(71) one finds that \ip Al ) <-> (-!)"" Wb x ), with 0£ = 1, By = 0; so \tp Al ) has vanishing 
amplitude on a sublattice and IV'bj) on the other. The two-body state \i>' A ^) cM>'b^) ~<y has occupation for spin 
a but not for spin —a on the site r = 0; under a lattice step translation it flips the spin and picks up a (-1) 
phase factor: \tp' Al ) ■<t\iPb 1 } -° IV'Si) 'flV'Ai) -o-; therefore it has double occupation nowhere and is a W = 
state (W = pair [38] [37]). 

The 4-dimensional irrep with basis ki = (n/2,n/2), k2 = (— vr/2, 7r/2), k3 = (n/2,—n/2) k.4 = 
(— 7r/2, — 7r/2) breaks into Ai © B2 © -E in C4,,; letting / = 1, 2, 3, 4 for the irreps A\, B2, E x , E y respectively, 
we can write down all the eigenvectors of (n r =o)ij = (k;|n r =o|kj), with i,j = 1, . . . ,4, as \tp'j) = Ojjk;}, 
where O' is the following 4x4 orthogonal matrix 



0' = 



1111 
1-1-1 1 
1-1 1-1 
-1-1 1 1 



(72) 



The state with non-vanishing eigenvalue is again the A\ eigenstate. After a little bit of algebra we have shown 
[57] that under Ti the subspace of A\ and B2 symmetry is exchanged with the one of E x and E y symmetry. 
Thus we can build a 4-body eigenstate of V^Hubbard with vanishing eigenvalue: \ip' Al ^>'b^i A^'e^'e^i -a ■ As 
before under a lattice step translation this state does not change its spatial distribution but a — » —a without 
any phase factor: W Ai ^'b 2 ) oWe x 4>e v }-<t <-> We x V'b„)<t|^a 1 4>'b 2 )-<t- 

Now we use these results to diagonalize n r =o on the whole set Shf (we could have done that directly 
by diagonalizing 6x6 matrices but we wanted to show the general method). The eigenstate of n r =o with 
nonvanishing eigenvalue always belongs to A\. The matrix n v has eigenvalues 3/8 and (5 times) 0, as predicted 
by Eq.(69). For r = the eigenvector of occupation 3/8 is |^ 0) ) = -^Wk x ) + The other Ai 

eigenstate of n T =o has eigenvalue and reads: |v4 } ~ y/^l^Ai) ~ ~j^\'4 , a 1 }- The other eigenvectors, whose 
symmetry differs from A u arc |^ 0) ) = |V>s 2 ), |^i 0) > = \i> Bl ), WT) = We x ) and |^ 0) ) = \ip' E ). One finds 



[57] that the transfer matrices T t of Eq.(71) such that \ip} e,> ) = J2j IVj ) T hn are 
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The reason why this choice of the basis set is clever is now apparent. The local basis at any site r splits into 
the subsets S a = {\ ( Pi > ' 1 }, 1^2°^)! I < y £> 3°' > ) } 5 an d St = {\tp^}, l^ ' 1 }, Iv^}}', a shift by a lattice step sends members 
of S a into linear combinations of the members of Sb, and conversely. 
Consider the 6-body eigenstate of K 



(74) 



In this state there is partial occupation of site r = with spin a, but no double occupation. It turns out that 
a shift by a lattice step produces the transformation 



\<S>af)<? < — > -]<S>af)- 



(75) 



that is, a lattice step is equivalent to a spin flip, a feature that we call antiferromagnetic property. Since the 
spin-flipped state is also free of double occupation, |$af)<t is a W = eigenstate. A ground state which is a 
single determinant is a quite unusual property for an interacting model like this. 
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Note that \ip[ ' ip\ ') is equivalent to IV'AiV'Ai); because this is just a unitary transformation of the A± wave 
functions; so j$AF)<T can also be written in terms of the old local orbitals (without any mix of the local states 
of different irreps of G): 

|*af) ct = \Tpl 1 1p'A 1 i>'B 2 }<T\lpB 1 ^E^E y }-<y (76) 

This form of the ground state lends itself to be generalised (see below). For Na > 4, k vectors arise that do not 
possess any special symmetry, the vectors -Rik are all different for all Ri G Cav, and we get eight-dimensional 
irreps of G. Recalling that \Shf \ = 2A?a — 2, one finds that Shf contains iV e = ^(-^r — 2) irreps of dimension 
8, one of dimension 4 and one of dimension 2 if Na/2 is even and N a = ^(-^p — 1) irreps of dimension 8 and 
one of dimension 2 if Na/2 is odd. 

To extend the theory to general TVa, we note that these k vectors, since Rik are all different, are a basis of 
the regular representation of C\ v . Thus, by the Burnside theorem, each of them breaks into Ai © A2 © -Bi © 
B2 © E © E; diagonalizing n r= o and the point group characters on the basis of the m-th eight-dimensional 
irrep of G one gets one-body states \ip^}, where I stands for the Ca v irrep label, I=Ai, A2, Bi, B2, E x , E y , 
E' x , E' y ; here we denote by E' the second occourrence of the irrep E. The ground state wave function in TL 
for the half filled case is a generalized version of Eq.(76). For even Na/2, we have [57] 

m— 1 m— 1 

with a =|, j. For odd Na/2, |$af)<t is similar but the maximum m is N and the |i/>'} states do not occour. 
\&af)ct is a W — state, transforms into — |<J?af)-<7 for each lattice step translation and manifestly shows 
an antiferromagnetic order (antiferromagnetic property). Since Wmibbard is a positive semidefinite operator 
|<&af) ct is actually a ground state. In the basis of local states these are the only two determinantal states 
(cr =|, J.) with the above properties. 

A few further remarks about |"J>af)<t are in order. (1) Introducing the projection operator Ps on the spin 
S subspace, one finds that P S |<I> af) ct = \<t>AF)<? / ,VS = 0,... , JV A - 1. Then, CT ("I>AF|H / Hubba r d| < I>AF) CT = 
X^s^i -1 ""(^afI W^Hubbard I "3? af)ct = 0, and this implies that there is at least one ground state of ^Hubbard in 
TL for each S. The actual ground state of -ffijubbard at weak coupling is the singlet |$af)<t- (2) The existence 
of this singlet W = ground state is a direct consequence of the Lieb theorem [58] . Indeed the maximum spin 
state I'&af -1 ) ' ' s trivially in the kernel of IFmibbard; since the ground state must be a singlet it should be 
an eigenvector of WHubbard with vanishing eigenvalue. (3) The above results and Lieb's theorem imply that 
higher order effects split the ground state multiplet of i^Hubbard and the singlet is lowest. 4) The Lieb theorem 
makes no assumptions concerning the lattice structure; adding the ingredient of the G symmetry we are able 
to explicitly display the wave function at weak coupling. 

Using the explicit form of Ps=o one finds that Ps=o| ( 1?af)ct = — Ps=o\®af)-<j- This identity allows us to 
study how the singlet component transforms under translations, reflections and rotations. In particular the 
antiferromagnetic property tells us that the total momentum is K to t = (0,0). To make contact with Ref. [59] 
we have also determined how it transforms under the d v operations with respect to the center of an arbitrary 
placquette. It turns out [57] that it is even under reflections and transforms as an s wave if Na/2 is even and 
as a d wave if Na/2 is odd. 

In the next Section we use these results, togheter with the non-perturbative canonical transformation, to 
study the doped 4x4 lattice at half filling. Since the non-interacting ground state at half filling is now well 
known and unambiguously defined the expansion (48) can be performed in a unique way. 



F. Pairing in the Doped Hubbard Antiferromagnet 

Let us represent the 4x4 lattice as in Figure 9 (A). 



5 6 7 8 
9 10 11 12 
13 14 15 16 



5 14 8 

6 2 3 7 
10 14 15 11 
9 13 16 12 



(A) 



(B) 



FIG. 9.: (A) The 4x4 lattice. The numbers in the squares enumerate the sites. (B) The 4x4 lattice after the 
dynamical symmetry operation. 



28 



Q 


Ci 


Cz 


Cz 


C 4 


C B 


C 6 


C 7 


Cs 


Cg 


C10 


Cn 


C12 


C13 


C14 


ClB 


C16 


Cn 


C18 


C19 


C20 


A 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


A 1 


1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


-1 


-1 


1 


1 


1 


1 


1 


-1 


-1 


B 2 


1 


1 


-1 


-1 


1 


1 


1 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


-1 


B 2 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


1 


-1 


1 


1 


-1 


-1 


-1 


1 


Ti 


2 


2 


-2 


-2 


2 


2 








2 














-2 


-1 


-1 


1 


1 








r 2 


2 


2 


2 


2 


2 


2 








2 














2 


-1 


-1 


-1 


-1 








Si 


3 


3 


3 


3 


3 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 














1 


1 


s 2 


3 


3 


3 


3 


3 


-1 


1 


1 


-1 


1 


1 


1 


1 


-1 














-1 


-1 


s 3 


3 


3 


-3 


-3 


3 


-1 


-1 


-1 


-1 


1 


1 


-1 


1 


1 














1 


-1 


s 4 


3 


3 


-3 
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3 


-1 


1 


1 


-1 


-1 


-1 


1 
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-1 


1 


Ai 


4 


-4 


-2 


2 








-2 


2 





-2 
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-1 


1 


-1 


1 








A 2 


4 


-4 


-2 


2 








2 


-2 





2 


-2 











-1 


1 


-1 


1 








A 3 


4 


-4 


2 


-2 








-2 


2 





2 


-2 











-1 


1 


1 


-1 








A 4 


4 


-4 


2 


-2 








2 


-2 





-2 
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-1 


1 


1 


-1 








«i 


6 


6 








-2 


-2 


-2 


-2 


2 
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6 
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-2 
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2 
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n 3 


6 


6 








-2 


2 








-2 


-2 


-2 
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6 


6 








-2 


2 








-2 


2 


2 





-2 























ni 


8 


-8 


-4 


4 
































1 


-1 


1 


-1 








n 2 


8 


-8 


4 


-4 
































1 


-1 


-1 


1 









TABLE X.: Character Table of the Optimal Group Q of the 4x4 model. 



We consider a 4x4 lattice one-band Hubbard Model and periodic boundary conditions. The one particle 
energy spectrum of the kinetic term K has 5 equally spaced levels having degeneracy 1, 4, 6, 4, 1 respectively. 
The space Group G containing the translations and the 8 C^ v operations is not enough to explain the degen- 
eracy 6. Indeed the largest dimension of the irreps of the space Group for this 4x4 lattice is 4. Nevertheless 
the 4x4 lattice is very special, since among the symmetry operations that preserve K there are some which 
are not contained in G and preserve ^Hubbard too. This fact implies that i/Hubbard is invariant under a new 
and largest symmetry Group; let us call it Q. As observed by previous authors [60] [61] the 4x4 lattice can be 
mapped in the 2x2x2x2 hypercube since each pair of next to nearest neighbor sites has two nearest neighbor 
sites in common. Due to the importance of the symmetry in our configuration interaction mechanism, below 
we explicitly build a new symmetry operations and derive the character table of Q [55]. 

1. Symmetry Group of the Ax A lattice 

Periodic boundary conditions are assumed and for example, the nearest neighbours of 1 are 2, 5, 4 and 13. 
Rotating the plaquettes 1, 2, 5, 6 and 11 ,12 ,15, 16 clockwise and the other two counterclockwise by 90 degrees 
we obtain the effect of the dynamical symmetry, that we call d, see Figure 9 (B). This transformation preserves 
nearest neighbours (and so, each order of neighbours) but is not an isometry, and for example the distance 
between 1 and 3 changes. Thus, this symmetry operation d is a new, dynamical symmetry. Including d and 
closing the multiplication table we obtain the Group Q with 384 elements in 20 classes (like G) . The complete 
Character Table of Q is shown in the Table X. As the notation suggests, the irreps A\ and Ai both reduce to Ai , 
in C\ v , while B2 and B2 both reduce to Bi- In the Table XI we report the one-body eigenvalues of for t = — 1, 
the degeneracy and the symmetry of each eigenvector. For pairs, the W = theorem of Section II A 2 ensures 
that no double occupancy is possible in the irreps Ai, B2, Ti, T2, Si, £2, S3, S 4 , A2, A3, Qi, SI2, IT 
and II2 not represented in the one-body spectrum. The exact interacting ground state of the system with 4 
holes is threefold degenerate and belongs to an irrep of Q that contains the irreps A\ and B\ of C 4 „ once and 
two times respectively [54]. This irrep is not contained in the one particle spectrum of K, which does not 
admit degeneracy 3. By the W = theorem, pairs belonging to it must be W = pairs. These W = pairs 
arise from the single-particle states of K with eigenvalue -2, and their irrep is contained in the square of the 
irrep of the single-particle level. 

The m states of the expansion (48) are all the W = pairs belonging to the ground state irrep. The binding 
energy |A(4)| can be computed analitically by truncating the configuration interaction expansion (48) to the 
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Energy 


Irrep of Q 


Degeneracy 


4 


B 2 


1 


2 


A 4 


4 





ft 4 


6 


-2 


Ai 


4 


-4 


Ax 


1 



TABLE XL: One-body spectrum for t = — 1. 

a states and also by means of exact diagonalization [54]. The results are listed below for different values of 
the on-site interaction U (energies are in eV) : 





U =0.1 


U =0.5 


U =0.7 


U =1 


A(4) cxact 


-0.053 


-0.32 


-0.67 


-1.18 


A(4) ana i ytic 


-0.078 


-1.95 


-3.83 


-7.81 



As expected, with increasing U the difference |A(4) cxact — A(4) ana i ytic | increases because the renormalization 
induced by virtual electron-hole exitations becomes important and it is no longer a good approximation to 
consider the a states only. Nevertheless the truncated canonical transformation still predicts the right sign of 
A. 

As observed in Section HIE the above canonical transformation applies when two holes are added to a 
non-degenerate vacuum. To study the system at and close to half filling, we have to use the results of Section 
III E. In the following we will solve the problem of two electrons added to the half filled system. 

2. W — Pairs and Quartets 

Here we use the antiferromagnetic ground state |$af°) to predict the possible symmetries of the doped half 
filled system. With 12 holes, in the U — » limit, there are two in Shf, the first-order ground states correspond 
to W = pairs. The symmetry of these W = states can be determined a priori from the Q irreps of Table 
XL Apart from the filled shells, two holes go to the fi 4 level. From the character of Q one can derive that 

Qt = Ai + B 2 + fi 4 + Li + L 2 + E 2 + E 3 + fti + n 2 + fi 3 ; (78) 

since the first 3 entries are present in Table XI, the W = theorem ensures that Li, r 2 , E 2 , E3, Qi, f2 2 and 
O3 pairs have no double occupation. It turns out that the spin and orbital symmetries are entangled, i.e. 
some of these pairs are triplet and the rest singlet. We can see that by projecting the determinantal state 
cj^CpJO) with k, p G Shf on the irreps not contained in the spectrum. One obtaines singlet W = pairs for 
Ti, r 2 , E 2 , fli. Other irreps yield W = triplet pairs. They are the three times degenerate irrep E3 and 
the two sixfold set of f2 2 and f2 3 symmetry. The above a priori argument hardly applies to the symmetries of 
W — quartets, because Of contains almost every symmetry and we do not know any W = 0-like theorem 
for quartets. However, we can still build the projection operators by Mathematica; we can project the 225 
quartets on the irreps of Q and carry on the analysis in an efficient, if not elegant, way [55]. It was found 
that the singlet W = quartets are 13 as many as the singlet W = pairs, and belong to the same irreps 
Fi, T 2 , E 2 , fii. Therefore, these are the possible symmetries of the first-order ground states with 14 holes. 
Exact diagonalization results [51] show that for U/t < 3 and 16 — 2 = 14 holes the ground state is sixfold 
degenerate, with a doublet of states with momentum (7r,7r) and a quartet with momentum (±7r/2, ±7r/2). 
It can be shown [55] that the computed ground state corresponds to an fii electron pair over the half filled 
system. For U/t > 3 and the same number of holes a level crossing takes place: the ground state is threefold 
degenerate and contains a state with momentum (0,0) and a doublet with momentum (ir,0) and (0, 7r). The 
computed ground state must be assigned to a E 2 electron pair over the half filled system [55]. In both cases, 
the symmetry of the ground state corresponds to a W = pair. 
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3. Pairing mechanism 



We consider the ground state of the 4x4 model with 14 holes; aside from the 10 holes in the inner A\ and 
Ai shells (see Table XI) the outer Q4 shell contains 4 holes in a W — quartet. We are in position to show 
that pairing between two electrons added to the antiferromagnetic 16-holes ground state (half filling) comes 
out. 

We recall from Section IIIF2 that, by comparing with exact diagonalization results [51], the ground state 
is assigned to fii at weak coupling and to £2 at a stronger coupling. We must preliminarily verify that Group 
Theory does not forbid obtaining these symmetries by creating W = electron pairs from the antiferromagnetic 
state |$af°}- This is the same as annihilating hole pairs. Since the state l^li? } i s a total symmetric singlet 
with vanishing total momentum, the labels of the quartets will be the same of the annihilated hole pairs. This 
operation can be done by hand, or with the help of Mathematica, and the answer to the preliminary question 
is adfirmative for fii and £2, but not for all pairs. One obtains 24 W = quartets out of the 28 states since 
the annihilation of I\ and F2 W — pairs gives identically zero. 

Next, we have to use the antiferromagnetic ground state I^aJ } as the non-interacting ground state of 
the configuration interaction expansion (48). In small clusters like the 4x4 one the one-body states are 
widely separated and the intra-shell interaction is much more important than the inter-shell one; therefore, we 
consider only the m states made removing two holes in Shf from l^&li? ), neglecting the high-lying unoccupied 
orbitals [32]. The explicit calculations [55] shown that the effective interaction is attractive for both fii and 
£2 W — pairs leading to a bound state with binding energy Aq 1 = —61.9 meV and Ae 2 = —60.7 meV for 
U = — t = 1 eV. Therefore, the true weak coupling ground state is a fii W = pair over the antiferromagnetic 
state. This result agree with exact diagonalization data [51]. 

IV. CONCLUSIONS 

We have presented the following evidence that the W = pairs are the quasi-particles that, once dressed, 
play the role of Cooper pairs: 1) as two-body states they do not feel the large on-site repulsion, that would come 
in first-order perturbation in any theory of pairing with any other kind of pairs. 2) The indirect interaction 
with the background particles gives attraction, and bound states with physically appealing binding energies. 
3) The same results are also borne out by exact diagonalisation in finite clusters, if and only if they have the 
correct symmetry and filling to give raise to W = pairs. 4) Both in clusters and in the plane, superconducting 
flux quantisation results from the symmetry properties of W — pairs. 

The setup of our theory of the effective interaction W e s between two holes (or electrons) is quite general; 
although we developed it in detail for Hubbard Models, we can readily include phonons and other ingredients. 
In principle we can obtain W e s by our canonical transformation, including systematically any kind of virtual 
intermediate states. The closed-form analytic expression of W c e we obtained includes 4-body virtual states. 
This describes repeated exchange of an electron-hole pair to all orders. One can envisage the pairing mechanism 
by spin-flip exchange diagrams that are enhanced by the Cav symmetry. The argument does not depend in any 
way on perturbation theory, and the equations retain their form, with renormalized parameters, at all orders. 
The previous exact-diagonalisation results of cluster calculations are special cases. The resulting integral 
equation (64), with the effective interaction (63), is valid for the full plane. Since an analytic treatment is 
prohibitive, we resort to a numerical treatment. We find that in the three-band Hubbard Model, 1 A2 pairs 
are more tightly bound close to half filling, but 1 E>2 pairs are favored when the filling increases. We remind 
here that these symmetry labels are not absolute, but depend on the choice of a gauge convention. We get 
attraction and pairing at all fillings we have considered (above half filling), but the binding energy of the 1 A2 
pairs drops by orders of magnitude as the filling increases; thus, there is no chance of superconducting flux 
quantization too far from half filling. So, we do not predict that superconductivity occurs outside some range 
of hole concentration. However, pairing is still there, even for large doping. The three-band Hubbard Model 
might be too idealized to allow a detailed comparison with experiments; however we stress that the approach 
presented is far more general than the model we are using, and can be applied to more realistic Hamiltonians. 

The above non-perturbative canonical transformation allows to study pairing fluctuations when the non- 
interacting ground state is unique (non-degenerate). Nevertheless, at exactly half filling, the ground state is 
highly degenerate and the configuration interaction expansion of Section IIIC is not well defined any more. 
To overcome the problem, we have removed the degeneracy in first order and with the help of the Lieb's 
theorem we have found the exact ground state of the half filled one-band Hubbard Model at weak coupling. 
We have shown that its symmetry quantum numbers are the same as in the strong coupling case. This result 
was estabilished for square lattices of arbitrary size. 
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Once the degeneracy is removed the non-perturbative canonical transformation can be applied and the 
interesting question of what happens by doping the Hubbard antiferromagnet can be explored analytically. 
Therefore, we test the W = pairing mechanism within the one-band Hubbard Model with periodic boundary 
conditions using exact numerical data on the 4x4 square lattice [51]. In particular, several workers have 
recognised that those data could be qualitatively understood by a weak coupling analysis [35] [36]; however 
the evidence of pairing was not seen and the symmetry analysis of the system was not complete. Furthermore, 
the role of W — states had not been discovered and the reason of the apparent success of the weak coupling 
analysis was not clarified. In our approach, the criteria exploited in the fully symmetric clusters of the Cu02 
plane are used to unambiguosly diagnose pairing. One more time the effective attractive interaction comes 
out from a spin flip diagram. This result lends further support to the general approach which predicts the 
existence of pairing fluctuations for systems of arbitrary size. 

On the other hand, the above results also prove that important ingredients are still missing and must 
be included. The 4x4 model shows evidence of bound pairs of non-vanishing momentum, in degenerate 
representations. This opens up the possibility of charge inhomogeneities and Jahn- Teller distorsions, that 
remain to be explored. 
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